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ABSTRACT 

We model gas and dust emission from regions 0.3—20 AU from a central low 
mass star in intermediate-aged (~10 7 years) disks whose dust is fairly optically 
thin to stellar radiation. The models treat thermal balance and chemistry self- 
consistently, and calculate the vertical density and temperature structure of the 
gas in a disk. The gas and dust temperatures are calculated separately. The 
models cover gas masses 10~ 3 — 1 Mj, dust masses 10~ 7 — 10~ 4 Mj, and treat 
solar type (G and K) stars. We focus on mid-infrared and far-infrared emission 
lines from various gas species such as the rotational lines of H2, OH, H2O and CO 
molecules and the fine structure lines of carbon, oxygen, sulfur, iron, and silicon 
atoms and ions. These lines and the dust continuum are observable by the Spitzer 
Space Telescope, and future missions including SOFIA and the Herschel Space 
Observatory. We find that the [SI]25.23/mi line is the strongest emission line for a 
wide range of disk and stellar parameters, followed by emission from [SiII]34.8/xm, 
[FeII]26/im, and [OI]63yum. [Fel]24/xm and rotational lines of OH and H 2 are 
strong when gas masses are high (> 0.1 Mj). Emission from the rotational lines 
of H 2 is more difficult to detect, unless disk gas masses are substantial (> 0.1 Mj). 
For emission from H 2 lines to be observable and yet the dust be optically thin in 
stellar light, the ratio of gas to small sub- millimeter sized dust particle mass in 
the disk needs to be > 1000, or at least an order of magnitude higher than that 
in the interstellar medium. This may be possible at intermediate stages in disk 
evolution, such as in the gas gathering stage of the core accretion scenario for 
giant planet formation, where most of the dust has coagulated into larger objects 
(3> 1mm) but the gas has not yet fully dispersed. Whereas the absolute fluxes 
observed in some lines such as [Fel]24/xm and H 2 S(0) 28/im primarily measure 
the gas mass in the disks, various line ratios diagnose the inner radius of the gas, 
and the radial temperature and surface density distribution of the gas. Many 
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lines become optically thick and/or suffer chemical or thermal effects such that 
the line luminosities do not increase appreciably with increasing gas mass. We 
predict it may be difficult for the Spitzer Space Telescope to detect < 1 Mj of 
gas in optically thin (in dust) disks at distances > 150 pc. The models presented 
here will be useful in future infrared studies of the timescale for the dispersion 
of gas in a planet-forming disk, and testing core accretion models of giant planet 
formation. 

Subject headings: infrared: spectra — planetary systems: protoplanetary disks 
— stars: circumstellar matter — stars: formation 



1. Introduction 

Circumstellar disks evolve from gas-rich structures with 1% of the total mass in dust to 
disks which are observed to be completely devoid of gas and dominated by emission from dust 
particles. During this process of evolution, the disk, which is initially very optically thick in 
dust with radial midplane optical depths in the visual ry > 10 3 , becomes optically thin with 
ry~l in a few million years (Haisch, Lada & Lada 2001). At this intermediate stage of disk 
evolution planetary formation is probably underway with dust particles colliding and growing 
to form larger objects (which greatly reduces ry), and with some residual gas. Disks with 
ages > 10 7 years are observed to have almost no gas (Zuckerman et al. 1995; Duvert et al. 
2000). Smaller dust particles in these disks have either dispersed or have coagulated to larger 
sizes, and these disks are very optically thin in dust with ry < 10~ 2 . These very optically 
thin disks are called "debris disks", as the dust in these disks is generally not primordial but 
is continually generated debris from colliding rocky bodies and planetesimals. 

The evolution of the dust component of disks has been extensively studied (e.g., Wei- 
denschilling 1977; Weidenschilling et al. 1997) in models where the dust grains in disks 
grow in size via collisional agglomeration. As dust grains grow and planetesimals eventually 
form, the disk becomes optically thin to stellar radiation and the fractional dust luminosity, 
fo = Lm/Lboi, the ratio of the disk infrared luminosity to the stellar bolometric luminosity, 
is observed to decrease with the age of the disk as t~ L76 (Habing et al. 1999; Spangler et al. 
2001; Zuckerman 2001). This decrease in dust luminosity occurs as micron-sized or smaller 
particles are lost through the effects of radiation pressure, Poynting-Robertson drag and 
collisions (e.g., Backman & Paresce 1993), and the population of larger bodies which are the 
ultimate source of the dust dwindles. 

Gas evolution in disks is less well understood. CO observational studies seem to indicate 
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that gas beyond about 30 AU disappears rather abruptly at disk ages > 10 7 years (e.g., 
Zuckerman et al. 1995; Duvert et al. 2000). Theoretical studies also indicate short gas 
dispersal timescales (Hollenbach, Yorke & Johnstone 2000; Clark, Gendrin & Sotomayor 
2001). However, there is some controversy regarding the presence or absence of gas in 
disks of ages ~ 10 7 years. Seemingly contradicting the CO observations, there may have 
been detection of ~0.1 — 1 Mj of molecular hydrogen gas in three younger (ages > 10 7 
years) debris disks by Thi et al. (2001) using ISO. These ISO results are controversial (see 
LeCavalier des Etangs et al. 2001; Richter et al. 2002; Sheret et al. 2003; Brandeker et 
al. 2004). If the ISO observations are confirmed, then perhaps CO is not a reliable tracer 
of gas mass in such systems, and may be subject to chemical effects such as freezing onto 
cold dust grain surfaces in the outer disks and optical depth effects in the inner disk. On 
theoretical grounds, although gas in the outer regions (r > 10 — 100 AU) may be subject 
to photoevaporation due to stellar radiation, gas inside these radii does not photoevaporate 
because it is more tightly bound gravitationally to the star. Viscous accretion onto the 
central star and viscous spreading to the outer photoevaporating disk likely dominates the 
gas dispersal of the inner (r < 10 — 100 AU) disk (Johnstone et al. 1998; Hollenbach et al. 
2000; Clarke et al. 2001). However, the timescale for viscous dispersal at r ~ 10 — 100 AU 
is uncertain. 

A determination of gas dispersal timescales in disks is important for understanding the 
formation of planetary systems. Two competing theories exist for the formation of gas giants. 
Gaseous Jovian-type planets are believed to form in disks via accretion of gas onto rocky 
cores of a few earth masses (e.g., Lissauer 1993; Pollack et al. 1996; Kornet, Bodenheimer 
& Rozyczka 2002) or by gravitational instability in disks leading to the formation of clumps 
which subsequently contract to form giant planets (e.g., Boss 2003). The presence or absence 
of gas in intermediate-aged or transitional disks (ry ~ 1, ages ~10 7 years) can potentially 
discriminate between these two theories. In the gravitational instability scenario, gas giants 
form quickly and the gas may dissipate rapidly. Longer gas disk lifetimes(> 10 6 years) 
facilitate the formation of gas giant planets through conventional core accretion models, 
allowing sufficient time for the building of a rock/ice core of a few M e and subsequent 
accretion of gas. At the gas accretion stage, there would be little dust (or equivalently, 
particles small enough to significantly radiate at infrared/sub-millimeter wavelengths and be 
observed) in the gas giant regions of disks, and substantial amounts (> 1 Mj, or a Jupiter 
mass) of gas. Core accretion theory suggests that systems may reside in this state for a 
significant length of time (few Myr, e.g., Pollack et al. 1996). At this stage, the ratio of gas 
mass to dust mass (where dust is defined as particles with sizes < 1 mm) would initially 
increase from the interstellar value of 100 before the final gas dispersal reduces the ratio to 
values <C 100 characteristic of the older debris disks. The detection of gas infrared emission 
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during this epoch of gas accretion onto rock/ice cores would help support the core accretion 
scenario. Conversely, if very little gas (<Cl Mj) is present at the moment when the dust 
disk first becomes optically thin (signalling the buildup of several earth mass cores), then 
the core accretion model is likely invalid. 

In the terrestrial planet region of the disk, the residual gas content of the disk at the 
epoch when embryos or protoplanets assemble to form terrestrial planets helps determine 
the ultimate mass and eccentricity of the planet, and therefore the eventual habitability 
of the planet (Agnor & Ward 2002; Kominami & Ida 2002; Kominami & Ida 2004). In 
particular, a finite range of gas masses produces earth-mass planets on circular orbits. If 
the gas mass inside 3 AU at ages of 10 65 to 10 7 years is 3> 10~ 2 Mj, the tidal interaction 
of the planets with the gas is sufficient to circularise the orbits of lunar-mass protoplanets, 
making it difficult for collisions between them to build earth-mass planets. On the other 
hand, if the gas inside 3 AU is <C 10~ 2 Mj, earth-mass planets can be produced, but in 
orbits considerably more eccentric than that of the Earth. 

At later debris disk stages, t 3> 10 7 years, the presence of even small amounts of gas 
(~10~ 2 Mj) can affect dust dynamics because of the drag forces they induce and, therefore, 
this gas may affect the disk structure significantly. This effect weakens the interpretations 
of gaps and ring signatures observed in dusty disks as being due to the presence of planets, 
as the presence of gas could give rise to the same structures (Klahr & Lin 2001; Takeuchi & 
Artymowicz 2001; Takeuchi & Lin 2002). 

In the work presented here, we model gas emission from a disk that has optically thin 
or marginally optically thick (ry < 5) dust, 1 so is perhaps at a stage where considerable 
rock/ice planetesimal and planet formation has already taken place, but may still retain 
appreciable amounts (from ~10~ 3 Mj to > 1 Mj) of gas left over from the rock/ice planet 
building phase. We predict expected emission from such a range of gas masses in optically 
thin dust disks as a function of the gas mass, dust mass, and stellar properties. The results 
will aid future observational programs that plan to look for gas emission in disks in order 
to understand the evolution of gas in disks, the timescales for gas dispersal, and the process 
of giant planet formation. New observational facilities such as the Spitzer Space Telescope, 
SOFIA and the Herschel Space Observatory are sensitive enough to detect infrared and sub- 
millimeter emission from small amounts of gas using tracers other than CO. These IR tracers 
may be better suited than CO for detecting gas in the 1 AU — 20 AU planet forming region 
around solar-type stars. The possibility of detection of gas emission from the planet-forming 



1 Henceforth, we shall call our models "optically thin" even though as discussed in §2.1, they apply for 
T v ^ 5, where ry is the radial optical depth through the midplane. 
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regions around intermediate-aged and evolved debris disks by future IR observational studies 
is one of the main motivations of the present paper. 

There has been extensive work on disk modeling of both young optically thick disks and 
older optically thin disks. However, most of this work has concentrated on either solving 
the disk structure by setting the gas and dust temperatures to be equal (e.g. D'Alessio 
et al. 1998, 1999; Dullemond 2002; Dullemond, van Zadelhoff & Natta 2002; Lauchame, 
Malbet & Monin 2003), or has assumed a gas density and temperature distribution and 
solved for the chemistry (Willacy et al. 1998; Aikawa et al. 1999, 2002; Markwick et 
al. 2002; Ilgner et al. 2004). We solve for both the gas and dust temperatures and the 
chemistry simultaneously in our models, and use the gas temperature to calculate the vertical 
density profile. A similar theoretical study has been made by Kamp & van Zadelhoff (2001), 
where the authors attempted to explain the observed H 2 emission from disks around two A 
stars, Vega and (3 Pictoris as seen by ISO and the lack of CO emission from these disks. 
These authors calculate the gas temperature by balancing the different heating and cooling 
mechanisms, but do not solve for the vertical structure of the disk. Our work adds the 
computation of vertical structure and is complementary to their model, since we focus on 
the inner regions of the disk, within about 20 AU from the star where the gas and dust are 
both warm (~100 K), whereas Kamp & van Zadelhoff (2001) considered disk regions from 
~40 — 200 AU. The inner regions are more interesting from the planet formation point of 
view, as many planets are expected to form close to the central star, within a few tens of 
AU if our solar system is typical. Furthermore, gas in regions far beyond this radius may be 
subject to photoevaporation from the central star and may be lost from the system early in 
its evolution. Our models are also complementary to the work of Glassgold, Najita & Igea 
(1997) and Glassgold, Najita & Igea (2004), who model the very innermost regions (r < 1 
AU) of young optically thick disks with emphasis on the problem of heating the surface 
layers and the near-infrared spectra. 

The paper is organized as follows. In the next section, we briefly describe some of the 
physical processes we consider. In §3, we describe our model in detail. We present the results 
of a fiducial case (§4) and then in §5 present results where we vary most of the parameters 
to see how disk emission is affected by their values. This is followed by a discussion of our 
results and conclusions (§6). 

2. Physical Processes in Disks 

Our disk models include chemistry and thermal balance in a self-consistent manner and 
calculate the density and temperature structure of the disk in both the radial and vertical 
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directions. The vertical structure of the disk is calculated by balancing the hydrostatic 
thermal pressure gradient with the vertical component of gravity from the central star. The 
inner (r^) and outer (r G ) boundaries of the disk and the distribution of the gas and dust 
surface density S(r) within these boundaries are input parameters. Our focus is on the IR 
emission detectable by the Spitzer Space Telescope, SOFIA or Herschel, and we therefore 
constrain our models to radial distances ~0.3 — 20 AU from the central star, where gas 
temperatures range from ~1000 K to 100 K respectively. The dust grains in this region 
are sufficiently warm to prevent ice formation on their surfaces, and therefore significant 
depletion of relatively volatile species such as gas phase O, C or S does not occur. We 
assume no significant opacity inside the inner edge at r\, so that there is an "inner hole". 

Although we focus on intermediate-aged and older debris disks which are relatively 
optically thin in dust opacity, the gas opacity can in fact be quite large at line center of 
strong transitions and even for some FUV continuum bands, which lead to photodissociation 
and photoionization of abundant gas species. It is instructive to calculate a fiducial radial 
(hydrogen) gas column density Nh in the equatorial plane. Assuming a gas surface density 
distribution E(r) oc r~ 3 / 2 , a vertical scale height of H = O.lr, a gas mass of ~0.1 Mj, and 
an outer radius r Q = 20 AU for a vertically isothermal disk, we find 



Since cross sections in the UV and X-ray for trace species can approach 10~ 17 cm~ 2 (or 
< 10~ 21 per H nucleus), we see that even small masses of gas can become optically thick in 
relevant bands. We also note that most of the opacity for radially declining surface densities 
is produced at r^. If stellar photons make it through the shielding at r^, they will generally 
penetrate the rest of the disk. 

We emphasize that our disk model is hydrostatic and our thermal and chemical solu- 
tions are steady state. Therefore, we implicitly assume that the timescales to reach chem- 
ical/thermal steady state are short compared to dynamic timescales, i.e., the timescale for 
radial or vertical migration of gas. Because of the high densities and subsequent rapid cool- 
ing timescales, the assumption of thermal balance is well justified. Chemical timescales, 
however, can be quite long in cool (T < 300 K) gas shielded from photodissociation. In such 
gas, cosmic ray ionization sets the equilibrium timescales to ~10 6 years, which is comparable 
to radial drift timescales. However, this work represents an important first step in predicting 
the chemical and thermal structure of the warm T ~100 — 1000 K, r ~ 1 — 20 AU regions 
of evolved disks and predicting their spectrum. 

We briefly describe below the chemical and physical processes considered in our models 
and refer the reader to the appendices for greater detail. 




(1) 
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2.1. Dust Physics 

Coagulation processes in intermediate-aged disks may lead to a distribution of solid 
particles which range in size from 1 /zm to 10 4 km (planets). We define "dust" to consist of all 
particles < a max = 1mm in radius. We assume that the size distribution of particles between 
a min and a max follows a power law n(a) oc a~ s , where 3 < s < 4 such that most of the surface 
area is in the smallest particles and most of the mass is in the largest particles. Observational 
studies of disks estimate the "dust" mass from emission at wavelengths ~ 0.4 — 2 mm. If the 
size distribution has the above power-law form, these observational studies are only sensitive 
to the mass of dust in particles with sizes < 1 mm. 

Implicit in these assumptions is that most of the mass of solids will be "hidden" in large 
particles, rocks, boulders and planetesimals which emit insignificant amounts of infrared and 
millimeter wavelength radiation. Thus, at intermediate disk ages (~ 10 7 years), where dust 
has coagulated up to planetary sizes but where significant quantities of gas may not yet 
have been dispersed, the gas to "dust" (meaning a < 1mm) mass ratio may increase above 
the value ~100 characteristic of interstellar conditions. For example, if 90% of the mass of 
the primordial dust in a disk grows by runaway accretion to cores/embryos of sizes 10 3 to 
10 4 km, and the remaining 10% of the mass is in solids with a distribution characterized 
by a min = 1/im, a max = 100km, and s = 3.5, the gas to "dust" mass ratio increases to 10 7 
without any gas or solid dispersion. Later, as the gas disperses, the gas to dust mass ratio 
may fall below the interstellar value. 

In order to express the large range of gas to (small) dust mass ratios and the absolute 
values of the gas or dust mass which may exist as intermediate-aged or older disks evolve, 
we consider a range of gas masses M gas ~ 10~ 3 — 1 Mj and a range of dust masses M dust ~ 
10~ 7 — 10~ 4 Mj (or fu ~ 10~ 4 — 10 _1 ). M gas and Mdust are treated as independent variables, 
so that we cover gas-to-dust mass ratios of 10 7 to 10. We emphasize that these masses lie 
at r < 20 AU. We do not consider gas masses smaller than 10~ 3 Mj because we find such 
small masses to be undetectable by the Spitzer Space Telescope, SOFIA or Herschel. We 
do not consider dust masses higher than 10~ 4 Mj because we find that the dust becomes 
significantly optically thick in the equatorial plane, and our model is valid only for small to 
moderate dust optical depths. 

Our models consider gas and dust separately but as a coupled system, and they are 
capable of treating the dust grain density distribution independently of the gas density 
distribution as a function of radius r and vertical height z. In addition, the grain size 
distribution can be arbitrarily varied with r and z. However, for simplicity and to limit the 
number of models to be presented, we make the following assumptions. 
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1. We assume that the dust mass density distribution p d (r, z) is proportional to the gas 
mass density distribution p(r,z). That is, for a given disk we assume that the same 
gas-to-dust mass ratio applies everywhere in the disk. This assumption means that 
we ignore settling by grains to the midplane. Turbulence may stir the gas and dust 
sufficiently to justify this assumption (e.g. Balbus & Hawley 1991). However, our main 
motivation is simplicity in the models; presumably by using an average (vertical) mass 
ratio of gas to dust, we will obtain a first order solution for the more complicated case 
of vertical differentiation. 

2. We assume that the grain size distribution does not vary with r or z. We consider 
our minimum grain size to be much larger than that in the interstellar medium. In 
intermediate-aged disks, the population of very small grains dwindles as coagulation 
processes lead to grain growth (e.g., Beckwith & Sargent 1996; Throop et al. 2001; 
Natta et al. 2004). In more evolved debris disks, small grains are continually replen- 
ished due to grinding down of larger objects by collisions, but radiation pressure forces 
act on sub-micron sized dust grains to rapidly remove them from the disk (Backman 
& Paresce 1993) on short timescales. Observations nevertheless seem to indicate the 
presence of small sub-micron sized particles (Li and Greenberg 1998) in a debris disk. 
We adopt a minimum grain size of a m i n = l//m for our standard case, but also consider 
other cases where we vary a min . In this paper, we only consider n(a) oc a -3 - 5 , i.e., 
s = 3.5. 

For a given mass in grains, the surface area oc a~^ 2 . Gas-grain collisional heating 
can be important in disks and due to their large surface area, small grains dominate the 
heating of the gas. Therefore, lowering a m i n raises the gas-grain coupling and usually raises 
the gas temperature in the disk. Lowering a m i n tends to raise the gas temperature for two 
other reasons. Smaller grains have higher efficiencies of heating the gas through the grain 
photoelectric heating mechanism (see Watson 1972; Weingartner & Draine 2001), which 
may dominate the heating once a m i n < 10~ 2 /im. More importantly for this paper, smaller 
grains are hotter than larger grains due to the heating by stellar radiation and the inefficient 
radiative cooling of small grains. Grain absorption and emission efficiencies are dependent 
primarily on their physical size relative to the wavelength of emission or absorption, and 
their composition. For our models, we choose silicate (amorphous olivine) grains and the 
absorption coefficients are calculated from a Mie code (SWolf, private communication) based 
on optical constants from the Jena database (Henning et al. 1999). We divide the grain size 
distribution into 32 logarithmically spaced bins in size a, and each size grain has a different 
temperature at a fixed point in space. 

We conclude this subsection with a discussion of the dust optical depths in the model 
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disks. Using the absorption coefficients discussed above, we calculate the optical depth Ty 
in the V band (0.5 /xm) of dust grains radially outward in the equatorial plane of the disk. 
For most of our models, Ty < 1. We do, however, run a few cases where the optical depth 
in the midplane Ty ~ 5. These disks are still optically "thin" in comparison to the younger, 
primordial dust disks where Ty ~ 10 3 — 10 4 . In our calculation of dust temperatures we ignore 
the heating of dust grains due to the infrared radiation field emitted by the surrounding dust, 
but we include extinction of starlight by dust grains along the line of sight to a particular 
point in the disk. Even at optical depths t v = 5, the dust grains are still primarily heated 
by attenuated starlight rather than by the re-radiated IR continuum field, so that our error 
in calculating the dust temperature is still small (< 10%), and our models are still valid. 
Typically, debris disks (t > 10 7 years) have Ty < 10~ 2 . However, intermediate-aged disks 
(t ~ 10 6 - 5 years) may have 10 -2 < t v < 5. 



2.2. Gas Chemistry 

The chemical composition of the gas disk is critical to the structure and evolution of the 
disk. Although H 2 cooling is often important in the radial range we consider, trace species 
such as S, O, CO and H 2 are often significant or even dominant coolants and for an accurate 
determination of the strength of emission lines, the abundances of atoms and molecules need 
to be calculated. Our treatment of disk chemistry is an adaptation of the Photodissociation 
Region (PDR) models of Tielens and Hollenbach (1985, hereafter TH85) and subsequent 
work (Kaufman et al. 1999). Elemental gas-phase abundances are taken from Savage & 
Sembach (1996) (Table 5, Appendix A). We use interstellar depletions of these elements, 
allowing for refractory dust but no ices in the r < 20AU disk region, where the dust is 
sufficiently warm to prevent ice formation. We have added sulfur chemistry and reactions of 
CH 3 and CH 4 , and there are a total of 73 species with atoms and ions of H, He, C, O, S, Fe, 
Mg, Si and molecules and molecular ions of H, C, O, S and Si (Table 6, Appendix A). To the 
PDR network of ion-neutral, neutral-neutral, and photoionization reactions, we have added 
X-ray photoreactions, and some high temperature reactions of importance under typical disk 
conditions. We have 537 reactions in all, with reaction rates in most cases taken from the 
UMIST Astro chemistry database (Le Teuff et al. 2000). A list of the added reactions, and 
the rates when they differ from the UMIST data, can be found in Table 7 (Appendix A). 

The formation mechanism of H 2 in disks may differ from the processes at work in the 
interstellar medium. In the ISM, H 2 forms mainly on the surfaces of cold dust grains, with 
Tdust < 25 K, typically. Dust in disks at 0.3 — 20 AU typically have temperatures T dust > 100 
K. It is not clear if H atoms would stick to such warm dust grains for sufficient time to allow 
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formation of H 2 molecules (Hollenbach & Salpeter 1971; Cazaux & Tielens 2002). We have 
conservatively omitted grain formation of H 2 . We find that other processes such as three- 
body reactions in dense regions and formation through H~ can completely convert atomic 
hydrogen to molecular hydrogen in self-shielded regions. We have also run cases including 
grain formation of H 2 and find that the resultant infrared spectrum is not significantly 
altered. 

Because the dust is optically thin, stellar and interstellar UV radiation is mainly at- 
tenuated in the disk by different gas species in various bands and lines. Large gas columns 
in the inner disk shield the outer regions from UV photodissociation caused by the central 
star. Similarly, the gas columns in the upper layers shield the midplane from the interstel- 
lar radiation field (ISRF). The self-shielding of H 2 and CO, which involves absorption of 
specific resonance lines rather than continuum radiation, is treated analytically as in TH85. 
Attenuation of continuum FUV radiation by atoms and molecules is determined by splitting 
the FUV stellar radiation field into nine energy intervals, and the photon flux in each inter- 
val is calculated from the stellar (and ISRF) spectrum. The nine intervals from 0.74eV to 
13.6eV include 0.74-2.6eV, 2.6-3.5eV, 3.5-4.3eV, 4.3-5. 12eV, 5.12-7.5eV, 7.5-10.36eV, 
10.36— 11. 26eV and 11.26— 13. 6eV where the lower and upper bounds in each case were cho- 
sen to correspond to dominant photoabsorption thresholds of species in each energy bin. 
The attenuated photon flux in each bin at a given position in the disk is dependent on the 
column densities of all the absorbing species towards the line-of-sight to the source. Pho- 
todissociation and photoionization reaction rates are thus accurately computed by taking 
into account the attenuation of photons in the energy range where absorption cross-sections 
peak. 

We also consider the effects of X-ray radiation from the star on the disk. X-ray luminosi- 
ties of ~10 -4 L bo i have been observed in stars as old as 10 7 — 10 8 years, and therefore include 
not only intermediate-aged disks but also a number of observed "debris disks" . X-rays affect 
the chemical structure of the disk through ionizations and heat the gas via collisions with 
secondary electrons generated by the X-ray produced primary electron. X-ray heating and 
chemistry may dominate in the inner regions of the disk and the surface layers where atten- 
uation by gas is low. Stellar X-ray spectra from weak-line T Tauri stars peak at about 1—2 
keV (e.g., Feigelson and Montmerle 1999) and the flux drops rapidly as the X-ray photon 
energy increases. Soft X-rays, with energies < 1 keV, do not penetrate significantly through 
a massive disk, their typical attenuation columns being about 10 21-22 cm -2 of gas. The 
penetration columns are larger by about two orders of magnitude for harder (3—10 keV) 
X-rays, but the stellar flux and photoabsorption cross-sections decrease sharply to make 
their influence on disk chemistry relatively unimportant in the shielded regions of the disk. 
On the other hand, in disks with low masses or lower column densities towards the line of 
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sight to the star (disks with large inner radii, for example), X-rays incident on the disk heat, 
ionize and dissociate significant numbers of atoms and molecules in the gas and our chemical 
code treats this chemistry in detail. For simplicity, we do not consider double or higher ion- 
izations, as in Maloney et al. (1996). The total X-ray photoabsorption cross section of the 
gas is taken from Wilms, Allen & McCray (2001), adapted to our gas-phase abundances (see 
Appendix B). We use the photoionization cross sections as originally published by Verner 
et al. (1993) and later updated by Verner & Yakovlev (1995). The primary ionization rates 
for each species are calculated for each species by integrating the product of the attenuated 
photon flux and the absorption cross section over the X-ray energy range (Maloney et al. 
1996). We also include secondary ionization rates and X-ray induced FUV radiation and 
the resulting photochemistry (Maloney et al. 1996). The details of our treatment of X-ray 
chemistry and heating are given in Appendix B. 



2.3. Heating and Cooling Processes 

One of the major innovations of this paper in terms of disk modeling is the separate 
calculation of gas and dust temperatures. Most previous models (with the exception of 
Kamp & van Zadelhoff 2001) assume that the gas temperature equals the dust temperature. 
While this is certainly true in the midplanes or central regions of very optically thick disks, 
it is not as true for debris disks, intermediate-aged disks, or the upper regions of younger, 
optically thick disks. Here, the relatively weak thermal coupling between gas and dust, 
and the different processes which heat and cool the gas and dust, can lead to significant 
temperature differences. This point was also made by Chiang & Goldreich (1997), who note 
that the gas temperature in the superheated layers of passive, optically thick dust disks can 
be significantly lower than the dust temperature due to molecular cooling. 

Gas in disks can be heated through many different mechanisms, and we consider the 
physical and chemical processes that are likely to dominate in these environments. The heat- 
ing of dust grains by the stellar radiation field and subsequent collisions with gas molecules 
transfers kinetic energy to the gas thereby heating it. Gas-dust collisions are the dominant 
heating mechanism for the gas through most of the disk, except in regions where the gas 
is heated to temperatures higher than the dust by other processes and here collisions with 
dust cool the gas. In regions where the columns through the disk are small enough to allow 
the penetration of FUV radiation, and H 2 is photodissociated to form atomic hydrogen, the 
heat associated with the formation of H2 raises the temperature of the gas. X-rays from 
the central star directly heat the gas and dominate the gas heating at the surface and in 
the inner disk. In addition to these mechanisms, we also include heating due to collisional 
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de-excitation of vibrationally-excited H 2 molecules, grain photoelectric heating, drift heating 
due to dust, exothermic chemical processes, photoionization of neutral carbon, and cosmic 
rays. All the heating mechanisms considered are described in some detail in Appendix C. 

The gas in the disk is mainly cooled through radiative transitions of the different species 
and, where the gas is warmer than the dust, through gas-grain collisions. We include cooling 
due to the fine structure and metastable lines of CI, CII, 01, Oil, Fel, Fell, SI, SII, Sil, Sill, 
the rotational lines of H 2 , H 2 0, OH, CO, CH, HCO + , and by Lyman a. We also include 
cooling due to the collisional ionization of atomic species, and vibrational cooling due to 
H2. Vibrational cooling due to CO was found to be unimportant for the disk parameters 
considered in this paper. For all the molecular rotational levels, except for H2 and CO, we use 
the analytical formulae of Hollenbach & McKee (1979, hereafter HM79). We derive new fits 
to their formulae for H 2 (ortho and para states) and OH (see Appendix D) using collisional 
rate data for 167 and 45 levels respectively. It was computationally expensive to include a 
detailed multi-level calculation for these molecules in the code, and we estimate the strengths 
of individual H 2 and OH lines using the derived density and temperature structure and 
solve for the level populations after a disk solution for the temperature, density and chemical 
abundance structure (i.e., the r,z dependence) has been obtained with the new analytical 
approximations for the total cooling by OH and H 2 molecules. We note that cooling 
transitions are often optically thick, and we use the escape probability formalism described 
in HM79. 

We solve for the rotational cooling due to H 2 in detail because one of our primary aims 
is to calculate the strength of the H 2 emission lines from the disk. We include the first 20 
levels of molecular hydrogen in our calculation of the cooling, and derive the population of 
each individually excited level. We have also added vibrational cooling of H 2 as in HM79 
and cooling due to collisional dissociation of H 2 molecules. We do a similar calculation of 
the first 20 rotational levels of CO and solve for the level populations to compute the total 
CO cooling. 

3. Disk Model 

In this section, we describe our disk model, discuss some of our assumptions and give 
the details of our numerical scheme for obtaining a disk solution. Because our main aim is 
to characterize observed emission from disks and infer disk properties from the comparison 
of observed spectra with our models, we have varied most of our input parameters across a 
wide range. However, we define a "standard case" where we assign typical values to each 
parameter. The results from the standard model are described in §4. 
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Surface density distribution We assume a power law distribution of surface density of 
both gas and dust with radius, S(r) oc r~ a , where r is the distance from the star. The 
nominal value for a is often taken to be 3/2 in the literature (Lynden-Bell & Pringle 1974; 
Chiang & Goldreich 1997). For our standard case, we assume a flatter surface density 
distribution defined by a = 1, motivated by observations of disks (Dutrey et al. 1996, Li & 
Greenberg 1998). The main effect of changing the slope of the surface density distribution 
is to alter the relative distribution of mass in the inner and outer disk. This changes both 
the slope of the dust continuum emission at longer wavelengths and the relative strengths 
of the gas emission lines that preferentially originate either in the warmer inner disk or the 
colder outer disk region. We study the effects of a change in the surface density distribution 
by varying a. As long as a < 2, most of the disk mass is in the outer regions; and for a > 
(and assuming a vertical scale height H that increases with r), most of the radial opacity 
occurs in the inner regions. 

Radial extent Our model disk is defined to have both an inner and outer edge. It is not, 
however, certain whether an inner (gas) hole would exist in evolved disks, even if the dust is 
heavily depleted in the inner regions. Theoretical models of dust disks often assume an inner 
truncation radius determined by the dust sublimation temperature, interior to which there 
is no dust (e.g., Bell et al. 1997; Chiang & Goldreich 1997; D'Alessio et al. 1998; Dullemond 
et al. 2002). SED modeling of observed dust disk emission also argues for the presence of 
inner (dust) holes in many disk systems (e.g., Natta et al. 2001, Dullemond et al. 2001). In 
evolved disks, where planetary formation is likely to have taken place, an inner gas cavity 
may be generated by massive planets that prevent outer gas from accreting onto the star, 
whereas viscous accretion of gas interior to the massive planets may rapidly drain gas there. 

The inner disk radius r\ is an important parameter of our models because it determines 
the penetration of stellar radiation into the disk and contains the hottest gas and dust. 
Emission from this region, therefore, is a dominant contributor to the higher frequency 
continuum from the dust and the higher excitation lines of the gas. Small values of imply 
higher densities and column densities at the inner edge and radiation can penetrate only a 
short distance into the disk before being attenuated. If disks with the same gas mass have 
larger r iy there is less rapid attenuation of stellar photons. Consequently, photodissociation 
and photoionization by X-rays and FUV radiation extends further into the disk and affects 
more of the gas mass. Observations of intermediate-aged and debris disks indicate disk inner 
radii ranging from a few tenths of an AU to tens of AUs (e.g., Liu et al. 2004, Calvet et al. 
2002, Li & Greenberg 1998). For our standard case we assume an inner radius of r\ = 1 AU. 
However, we include several cases where we vary the inner radius to see how this parameter 
affects disk thermal structure and chemistry and the resultant spectrum. 
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We also define an outer disk radius r Q . Our fiducial disk is truncated at r Q = 20 AU. 
Unlike r iy the exact choice of r Q does not affect our results significantly because any gas 
beyond r Q will be too cold to significantly contribute to the IR lines we focus on in this 
study. The choice of 20 AU is motivated by the following four reasons, (i) Disks may 
lose much of the outer gas through the process of photoevaporation early in their evolution 
(Hollenbach et al. 1994; Adams et al. 2004). (ii) Our interest lies in the planet-forming 
regions of the disk, and from our understanding of the formation of the solar system and of 
the evolution of disks, planet formation appears to take place within a few tens of AU from 
a star, (Hi) Glassgold et al. (2004) present models of inner (< 1 AU) disks whereas Kamp 
& van Zadelhoff (2001) present models of outer (> 30 AU) disks. We, therefore, present this 
complementary study, which is especially relevant to the Spitzer Space Telescope, SOFIA and 
Herschel since the gas and dust temperatures span a range (~ 50 K - 500 K) which produce 
copious IR emission in the 4 fim — 200 /im wavelength range. Mid-infrared emission from 
gas species at temperatures less than ~ 100 K is likely to be insignificant, so the neglect 
of any gas beyond r Q does not affect our emission line calculations, (iv) Gas molecules, 
such as H 2 0, may freeze on grain surfaces when the dust grain temperature drops below 
~ 80 — 100 K. We therefore truncate the disk at a radius where ice mantle formation is 
likely to become important on grains, i.e., when dust grain temperatures fall below ~ 100 
K. To summarize points (Hi) and (iv), the outer edge need not correspond to a real cut-off 
in surface density. We are interested in the IR emission from warm gas in the r < 20 AU 
region. Material beyond this outer radius will be either too cold or frozen on grain surfaces 
to radiate appreciably as IR line emission from gaseous species. 

Stellar Spectrum The stellar spectrum is chosen to be a modified Kurucz model, and 
we consider stars of two spectral types, a G star at 6000 K and a cooler K star with a 
temperature of 4000 K. The central star at these disk ages is expected to have a significantly 
higher FUV luminosity than its main-sequence counterpart and we account for this increased 
luminosity in our stellar models. An accurate determination of the FUV flux is especially 
important in solving for the disk chemistry as the FUV initiates numerous photodissociation 
and photoionization reactions. We follow a prescription by Kamp & Sammar (2003) for 
our age-dependent FUV fluxes, where the FUV flux scales inversely with the age of the 
star (Ayres 1997). We assume an age of 10 7 years for the G star in our models. IUE and 
rocket experiment data are added to a Kurucz photosphere model to thus obtain a semi- 
empirical age-dependent spectrum for a G2 spectral type star. We use chromospheric models 
by Cohen(private communication) for the star e Eri as our input stellar model for the K star. 
Figure 1 shows the spectra we use, the solid line indicates the G star and the dashed line 
the K star. 
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For the stellar X-ray spectrum, we use a broken power-law fit to ROSAT observations 
(Feigelson & Montmerle 1999) for a wTTS, where the luminosity per keV peaks at about 
2 keV. The luminosity per keV is thus given by L(E) oc E~ 1 ' 75 ,E > 2 keV and L(E) oc E 
for E < 2 keV. This X-ray spectrum is normalized to the stellar X-ray luminosity between 
0.5—10 keV. For our fiducial case, the X-ray luminosity is assumed to be 10~ 4 times the 
bolometric luminosity of the star. 

Numerical model The disk density and temperature structure is calculated iteratively 
from the midplane to the surface starting at r = V{ and then working radially outward in the 
disk. We start with an estimate of the gas density at the point (r = r,, z — 0) and compute 
the dust and gas temperatures. 

The dust temperature is calculated from thermal balance between the absorption of 
stellar radiation by the dust and subsequent re-emission in the infrared. Heating of a grain 
due to absorption of infrared radiation coming from surrounding dust particles is ignored 
in our optically thin dust disk models. The dust temperature (at a given size a) needs 
to be iterated only because the dust vertical distribution is determined by the gas vertical 
distribution, and the dust temperature is a function of z because of the varying distance to 
the star. 

The gas temperature is determined by thermal balance between the heating and cooling 
processes described earlier, and this then determines the gas pressure. The condition of 
vertical hydrostatic pressure equilibrium, 

dP _ -GM*[im H n H (r,z) 

dz ~ ( r 2 + ^2)3/2 ^ ) 

determines the pressure as a function of disk height, where P = /c B (n(H)+n(H 2 )+n(He))T gas , 
/i is the mean molecular weight, m H is the mass of a hydrogen atom, n# is the hydrogen 
nucleus number density, and n(H), n(H 2 ) and n(He) are the number densities of atomic and 
molecular hydrogen and helium respectively. At each height z from the midplane, the density 
and temperature combination that gives the pressure from Eq.(2) is found. The vertical grid 
is adaptive in nature and determined by the condition of pressure gradient balancing gravity. 
The vertical increments in height are chosen to correspond to equal factors of decrease in 
thermal pressure (a factor of 0.8 in most cases) and for regions where the temperature slowly 
varies with z, this also corresponds to a logarithmically decreasing column density in each 
cell. This ensures that the escape probabilities for the cooling transitions increase slowly and 
smoothly as we get farther from the midplane, and enables accurate determination of the gas 
temperature. We ensure that the density contrast from surface to midplane is at least four 
orders of magnitude, and the number of cells and their spacing are calculated to meet all 
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these criteria. After the entire vertical grid has been solved, the gas density is integrated over 
z and compared to the required gas surface density, S(r). The midplane density is re-scaled 
and the vertical temperature and density structure is thus iterated until the surface density 
converges to the required value. 

Once we have the solution at we move outward to solve at the next radial grid point. 
The radial grid is logarithmically divided in r, except for the inner disk which is more finely 
sub-divided. The column density towards the star is dominated by the dense inner regions, 
and this is where the gas absorbs most of the stellar FUV and X-rays. This attenuation 
needs to be calculated accurately, and we have divided the inner disk into cells of slowly 
increasing column density going outwards in r, with the first radial cell (at midplane) having 
a hydrogen gas column density 10 20 cm -2 , which corresponds to a small attenuation of the 
X-ray and UV fluxes. 

4. Standard Disk Model Results 

In this section, we describe in some detail the results of our fiducial disk model ("the 
standard case"), the input parameters for which are listed in Table 1. The standard model 
has a gas mass M gas = 10" 2 Mj and dust mass 1CT 5 Mj. In the next section we present a 
full parameter study varying both the gas mass and the dust mass, the main parameters of 
this work. The other parameters we treat as "minor" parameters, and in the next section 
we vary each of these parameters in turn, holding all other parameters at the fiducial values, 
and study its effect on disk structure and emission. These minor parameters include the 
disk inner radius, the minimum grain size, the stellar X-ray flux and the surface density 
distribution. 

Our fiducial gas to dust mass ratio of 10 3 is higher than the interstellar value, but, as 
we have discussed, this is possible at an intermediate stage where dust has coagulated but 
gas has not yet dispersed. We have chosen this value for our fiducial case primarily because, 
as shall be shown below, this ratio allows detection of gas emission lines over the continuum 
for the smallest gas masses with the dust midplane optical depth still being < 1. 

Our standard model disk solution results in a flared disk with an inner puffed region 
where the gas is subject to the unattenuated, strong stellar radiation field and heated to 
high temperatures. Dullemond et al. (2002) find a similar effect in younger, optically thick 
(in dust) disks. Figure 2 shows the disk surface where the density is 10 -4 times the density 
in the midplane and contours where the vertical column densities to the surface are 10 21 and 
10 22 cm -2 . Stellar UV and X-ray fluxes are absorbed at these column densities, and hence 
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in these regions. The hot inner edge or rim has a somewhat lower density compared to gas 
just beyond the edge because it is hotter and more extended vertically. This region shields 
much of the rest of the disk from UV and X-ray photodissociation and photoionization by 
starlight. The density increases slightly with radius progressing from this inner edge as the 
temperature falls steeply, and the disk scale height decreases. Recall that for a vertically 
isothermal disk the scale height H of the disk is given by H oc (r 3 T(r)) 1 / 2 (e.g., Adams, Shu 
& Lada 1988, Shakura & Sunyaev 1973). Therefore for a radial temperature gradient steeper 
than 1/r 3 , the scale height decreases with radius. Further into the disk, the temperature 
falls less steeply and here H / r now increases with radius, giving rise to a flared disk solution. 
The density falls less steeply with height as compared to an vertically isothermal solution 
in most regions of the disk, because the gas temperature in general rises with z. We first 
describe the main heating and cooling mechanisms in the disk and the resulting temperature 
structure and then discuss the chemistry and emission from the disk. 



4.1. Heating, Cooling and Temperature Structure 

Figures 3, 4, 5 and 6 show the heating, cooling and the radial temperature and density 
structure in the disk at the midplane and the vertical temperature and density structure in 
the disk at two typical radii, of 2 and 10 AU. Figure 3 shows the main heating terms for 
our standard case in the midplane of the disk, as a function of disk radius. Heating due 
to formation of H2 is high where hydrogen is predominantly atomic, but quickly decreases 
as self-shielding turns the gas molecular, at hydrogen nuclei column densities of ~ io 21 ~ 22 
cm -2 . X-rays are a strong heating agent in the unshielded innermost regions of the disk, but 
get attenuated rapidly due to the high midplane radial column density N(r) ~ 10 24 cm~ 2 
(see Eq.(l)). The penetration column for 2 keV X-ray photons, where the X-ray spectrum 
peaks, is ~ l/a(E) ~ 5 x 10 22 cm -2 which implies a penetration depth of ~ 0.003 AU. 
Heating due to X-rays becomes unimportant in the midplane for r > r; into the disk. 
X-rays, however, can penetrate further at the surface of the disk, where the attenuation 
columns of gas towards the line-of-sight to the star are low. X-rays are the dominant heating 
mechanism at the surface. Gas-grain collisions are the strongest gas heating source through 
most of the equatorial regions of the disk even though the dust surface area per H nucleus is 
much smaller than the interstellar value(~10~ 3 times that in the ISM for the standard case). 
Note that in the innermost regions, the gas-grain collisions become a cooling mechanism, 
rather than a heating mechanism, as the gas becomes hotter than the dust grains due to 
X-ray and H2 formation heating. Heating due to cosmic rays, grain photoelectric heating 
and drift heating are negligible in this disk model. Kamp & van Zadelhoff(2001) find drift 
heating to be the dominant heating agent in their models of gas disks around A stars, but 
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we find that our more careful calculation of drift velocities yields values that are typically 
too low (~10 3 cm s -1 or less) for viscous heating to be important in our ~20 AU disks. 
Nevertheless, we include this mechanism as it may become important in other regions of 
parameter space. The dominant coolants in the midplane of the disk (Figure 4) are OH, 
H2O, SI, Sill, H2 and Fell in the inner regions and 01, CO and OH in the outer disk. Above 
the midplane SI and 01 are strong coolants. At radii larger than ~ 7 AU, the IR continuum 
field from the warmer dust excites the upper level of the [SI] line, and subsequent collisional 
de-excitation leads to a gas heating rather than cooling by this line, as seen by the sharp 
cutoff in the SI cooling in figure 4. 

The gas temperature in the midplane of the disk (Figure 5) decreases with radius ap- 
proximately as T gas oc r -0 - 8 . The plot also shows the dust temperatures of the smallest 
and largest grains. These temperatures are higher than that of the gas over most of the 
disk interior and therefore gas-grain collisions heat the gas here. In the outer regions of the 
disk midplane, the gas temperature falls well below the dust temperature, as the decreasing 
density leads to a decrease in collisional coupling between grains and dust which lowers the 
heating rate of the gas. Figure 5 also shows the number density of gas in the midplane which 
falls approximately as 1/r 2 . 

Figure 6 shows the variation of gas temperature with height at two different disk radii, 
2 and 10 AU, and the gas number density with height at 2 AU, with the dominant heating 
sources and coolants marked at various heights. In the equatorial region where gas-grain 
collisions dominate the heating, the density(dashed line in Figure 6) can be seen to decrease 
rapidly with z (recall that for isothermal disks, n oc e~ z2 l H2 , where H is the scale height), 
leading to a reduced gas-grain coupling and heating, and causing a drop in temperature with 
height. Several scale heights above the midplane, the attenuating column to the star drops 
to sufficiently low values that X-ray heating begins to dominate and the gas temperature 
rapidly rises. The main coolant at the surface is [OI]63/im, which is optically thick at these 
column densities. The column to the surface decreases with height, and the line becomes 
more optically thin and more efficient at cooling. This leads to a drop in temperature again 
in the upper disk atmosphere. Figure 6 also shows the temperature profile at a larger disk 
radius (10AU), where gas temperatures are much lower. The main heating sources remain 
gas-grain collisions at the midplane and X-rays at the disk surface, but the main coolants 
in these regions are CO and 01, as the temperatures are too low (< 100 K) to significantly 
excite the transitions of H 2 , SI, Sill or Fell, and the densities and OH and H 2 abundances 
are too low for OH and H 2 to dominate the cooling. 
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Chemistry Figure 7 shows the abundances of various species as a function of height above 
the disk at a radius of 2 AU. We discuss below the essential species in the hydrogen, OH/H 2 0, 
C+/C/CO, Si+/Si, Fe+/Fe, and S+/S chemistry. 

In the innermost disk, for r < 1.1 AU, hydrogen is predominantly atomic. H 2 is formed 
via reactions of H with H" and three body reactions, but is rapidly destroyed by reactions 
with O to form OH and by photodissociation. In regions high above the midplane, where 
there is less attenuation of stellar UV, photodissociation of H 2 dominates. Beyond 1.1 AU, 
self-shielding prevents UV photodissociation of H 2 . The destruction reaction with O has 
a high activation barrier at AE/k = 3200 K and this reaction rate rapidly falls as the 
temperature decreases radially outward. Consequently, hydrogen turns fully molecular in 
the dense midplane at r ~ 2AU. At the surface of the disk at r = 2 AU, lower attenuation 
columns lead to photodissociation of H 2 , and hydrogen is atomic as is shown in Figure 7. 

Carbon turns molecular in the midplane even at the inner edge of the disk where hy- 
drogen is still atomic, and is in the form of CO throughout the disk midplane. In the upper 
regions of the disk, CO is photodissociated into atomic carbon, which gets ionized even higher 
up in the disk to form C + (Figure 7). CO is formed via the reaction of C with OH in the 
inner disk where the temperatures are high and OH is more abundant. O reacts with CH to 
form CO in the outer disk. The main destruction routes of CO are photodissociation by the 
stellar UV field at the inner edge of the disk and in regions above the midplane everywhere 
in the disk. In the shielded interior of most of the disk, all the available carbon forms CO, 
as the only destruction mechanism for CO molecules is collisions with the very few He + ions 
that are formed by cosmic rays. 

The gas phase elemental oxygen abundance is a factor of 2.3 higher than the carbon 
abundance (Table 5, Appendix A) and all the oxygen not locked up in CO is in the form 
of 01 throughout the disk. The stellar photons are not energetic enough to ionize oxygen, 
and trace amounts of the atomic oxygen form 2 , OH and H 2 0. In the inner disk, where the 
temperature is high (T > 300 K), O reacts with H 2 to form OH, which then reacts again 
with H 2 to form water. OH also reacts with O to form 2 throughout the disk, which is 
destroyed by photodissociation. H 2 and OH are the dominant coolants at the high densities 
and temperatures in the innermost disk. Water is photodissociated through the disk more 
readily than OH as it has a lower photodissociation threshold (~4.3 eV). In disks that are 
more massive, larger attenuation columns shield H 2 from photodissociation and due to its 
high abundance it dominates the cooling, as will be discussed later. As the temperature 
drops further into the disk, the production of OH and H 2 drops, because the reactions of 
O and OH with H 2 have high activation barriers. The resulting decrease in the abundance 
of OH lowers its strength as a coolant further into the disk. The destruction of OH is 
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mainly via photodissociation which rapidly decreases with disk radius as the optical depth 
in the relevant FUV band increases, and this keeps the abundance of OH from dropping 
sharply. The main reaction for the production of OH and H 2 in the cooler outer regions is 
recombination of H 3 + , which is initiated by cosmic ray ionization of H 2 . 

Sulfur is predominantly atomic throughout most of the disk, except very near the mid- 
plane where it is partly molecular and at the disk surface where it is ionized to form SII. 
In more massive disks, which shield sulfur molecules from photodissociation, we find that 
most of the sulfur is SO2 at regions near the midplane. Silicon and iron are ionized through 
most of the disk, and are neutral only in denser disks where there is sufficient shielding from 
photoionization. Limited silicon chemistry is included in our models and we find that silicon 
is ionized through most of the disk. For more massive disks with M gas > 0.1 Mj, silicon is 
atomic in dense regions such as near the midplane, but does not turn molecular even for our 
most massive disks. 

SI, Sill, Fell and H 2 S(0) all have their first excitation level ~500 K above ground and 
are excited under similar conditions. The gas needs to be warm, > 100 K for emission from 
these lines to be strong. Note that sulfur is a much stronger coolant than H 2 though both 
species have similar excitation energies. At disk gas densities, these coolants are in LTE 
so that their cooling is proportional to their abundance times their spontaneous transition 
probability. Though the abundance of sulfur is 4.5 orders of magnitude lower than that 
of H 2 , its spontaneous transition probability is 4.7 x 10 7 that of H 2 S(0). Thus, at high 
densities, and assuming the lines are optically thin, [SI]25.23/im should dominate H 2 28 /xm 
by a factor ~1000. In fact, the [SI]25/im transition is optically thick, so that the [SI]25/im 
line is typically ~10 times the H 2 line strengths. Cooling from SI decreases as the gas 
temperature decreases, and in the outer disk the most important coolants are CO and 01 at 
the midplane and OI at the surface where CO is photodissociated. 

Disk emission and spectrum Figure 8 shows a model spectrum for this standard disk 
at a distance of 30 pc, with the dominant lines in the 5—40 /im band, superimposed on 
the dust continuum. (The resolution of Spitzer-IRS in the high-resolution mode (R = 600) 
has been assumed for computing the emission spectra.) The lines shown in the emission 
spectrum are strong enough for nearby sources (< 100 pc) to be potentially observable by 
the Spitzer Space Telescope, or future infrared telescopes, such as SOFIA or Herschel. The 
IRS instrument on the Spitzer Space Telescope has a 3a flux detection limit of ~10~ 18 W 
m~ 2 , with an integration time of ~500s. 

The [SI]25.23/im line stands out as a strong coolant in our disk spectrum and as an 
important diagnostic of the presence of gas at ~1 AU in disks. The possible detection of gas 



- 21 - 



disks in sulfur is one of our main results, as in many cases this line dominates the spectrum 
and is stronger than the H 2 lines. The [OI]63/im and [SI]56.0/xm lines may also be strong 
enough for detection by SOFIA. Other detectable lines in this figure are the H 2 S(l) and 
S(2) lines, the [FeII]26.0/im line, the [SiII]34.8/xm line, and water lines at 33/zm (o-H20, 
66,1 — > 5 5)0 ) and 36/im (o-H20, 6 5j2 — > 5 4 ,i). 

Figure 9 shows a similar spectrum for a disk with an inner radius of 10 AU, which is 
seen to have no strong OH or H 2 lines, but has detectable H 2 S(l), H 2 S(2), [SI]25.23//m, 
[FeII]26/im and [SiII]34.8//m lines. The effect on the gas of changing the inner disk radius is 
discussed in more detail in §5.2. 

Table 2 lists the observable mid-infrared and far infrared lines with their wavelengths 
and the telescope instruments that can observe them. Gas emission lines are potentially 
important diagnostic tools of physical and chemical conditions of gas in disks and could be 
useful in detecting even very small amounts of gas. However, most of these lines are optically 
thick with the exception of H 2 S(0) and would correspondingly diminish in strength with 
the inclination of the disk, being the strongest for a face-on orientation. Molecular hydrogen 
remains optically thin in the S(0) 28/im line and the line strength is not affected by the 
orientation angle of the disk. Moreover, it is less dependent on the chemical network and 
the cosmic abundances assumed and directly traces the gas in the disk. However, the H 2 
S(0) line is weaker than other gas emission lines and hence is more difficult to detect. For 
pure sensitivity to the presence of small amounts of gas at r ~ 1 — 2 AU, the [SI]25.23/im 
line is predicted to be the best probe. 



5. Parameter Survey and Emission Line Detectability 

5.1. Variation of Spectra with M_ gas and M dust . 

We have conducted a parameter survey to explore the probability of gas emission line 
detection above the dust continuum emission for a range of gas masses and dust masses 
and for a range of stellar and disk parameters (Table 3). We limit the range of parameters 
to those which might give us detectable infrared line fluxes, and those whose dust optical 
depths Ty < 5. Figure 10 shows the change in IR line luminosities due to variations in 
gas and dust masses. For a disk at a distance of 30 pc, the flux detection limit of the IRS 
(high resolution) instrument on Spitzer for a 3a 500-second detection corresponds to a line 
luminosity of ~3 x 10~ 8 L (or 10~ 18 W m~ 2 ). For comparison, SOFIA at these wavelengths 
(10-28 fim) has a minimum detectable 3a line flux of ~10~ 17 W mT 2 (3 x 10~ 7 L Q ), and 
Herschel ~7 x 10~ 18 W mT 2 (2 x 10~ 7 L ) at its shortest observing wavelengths (60-200 
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/xm). Instruments on both telescopes have resolving powers ~10 3-4 . Figure 11 presents 
these results in an alternate way, which readily shows the parameter space where lines are 
detectable. 

The S(l) and S(2) lines of H 2 arise from gas with temperatures > 150-200 K, whereas 
S(0) emission has a significant contribution from slightly cooler gas (T > 100 K). At very low 
gas masses (10 -3 Mj), the midplane column density through the disk is ~ 10 23 cm -2 , and 
hence low enough to allow penetration by stellar UV and X-ray photons through much of the 
inner disk and to photodissociate H 2 . The low amount of warm H 2 gas results in very low line 
strengths of all three H 2 lines, as can be seen in Figure 10. As the disk gas mass is increased 
to 10~ 2 Mj, the column density increases allowing shielding of H 2 , increasing its abundance 
in the inner, warm regions of the disk. The gas in the inner disk where the emission from 
the S(l) and S(2) lines peak, is heated by collisions with dust, H 2 formation heating and 
by X-rays. The increased amounts of warm H 2 gas (> 100 K) result in a sharp increase in 
the line strengths. As the gas mass is further increased to 0.1 Mj, the gas temperature at a 
fixed radial column density initially decreases because the H 2 cooling rises faster than the 
X-ray heating with density. Thus, there is less mass and surface area of warm gas. As the 
gas mass increases to ~ 1 Mj, the gas heating by gas-dust collisions takes over and raises 
the gas temperature and the S(l) and S(2) luminosities. The S(0) line emission peaks at 
regions slightly further out where dust dominates the heating and it monotonically increases 
with gas mass. At low gas masses, the observed H 2 S(0) emission directly correlates with 
the mass of warm gas in the disk, regardless of the dust or total gas masses, and the warm 
gas (T > 100 K) mass in a disk at 30 pc needs to be above ~10~ 3 Mj for detection by the 
Spitzer Space Telescope. At higher gas masses, ~1 Mj, dust collisions are the main heating 
mechanism that heat H 2 gas. In these relatively massive disks, H 2 line emission does not 
increase proportionately with mass, as the lines start to become optically thick (at disk gas 
masses > 0.05 Mj for S(l) and S(2), > 0.2 Mj for S(0)). Figure 10 shows that for line 
luminosities > 5 x 10~ 8 L , the luminosity generally scales with a low power(< 1) of gas 
mass. This means that for disks at distances d 3> 30 pc, considerably more gas mass is 
required for detection than if L oc M gas is assumed. This important point will be discussed 
more in §6. 

An increase in total dust mass of the disk produces increased emission due to an increase 
in heating of the gas via gas-dust collisions. Dust collisions are an important heating mech- 
anism in the bulk of the H 2 emitting regions. As the dust mass is decreased, the decreased 
coupling between gas and dust lowers the gas heating rate and the gas temperature and 
therefore the emission. On the other hand, dust produces IR continuum emission which can 
obscure the detection of infrared lines (typically, systematics introduce baseline variations 
which make it impossible to detect lines at the < few per cent continuum level). There is 
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thus an optimal gas-to-dust mass ratio ~ 10 3 — 10 5 and gas mass > 1CT 2 M j (for disks at 30 
pc) where H 2 emission can be detected. Given > 10~ 2 Mj of total gas mass, a gas to dust 
mass ratio < 100 produces a line to continuum ratio of < 5% for the highest Spitzer Space 
Telescope resolving power (R=600), rendering the line difficult to detect above the strong 
continuum. On the other hand, if the dust abundance is lowered to the point where the 
gas to dust ratio > 10 5 , then the grain-gas collisional heating is so drastically reduced that 
the gas is cold, and the H 2 lines are too weak to detect in the absence of other gas heating 
mechanisms. 

The SI, Fell, and Sill emission line strengths are complicated by chemistry in the disk. 
For disks of low gas mass (10~ 3 — 10~ 2 Mj), the main sulfur species in the denser regions 
near the midplane of the disk is atomic sulfur, and sulfur is photoionized near the surface. 
Iron and silicon are completely ionized throughout such disks. In this low gas mass limit, 
the emission increases as the gas mass increases and as the dust mass increases. A higher 
dust mass, in general, implies more heating through gas-grain collisions and hence warmer 
gas. As the gas mass increases to 0.1 Mj or higher, the main sulfur-bearing species near 
the midplane is SO2, as sulfur turns molecular at these high densities and opacities. Here 
most of the atomic sulfur is in a layer above the midplane where there is enough penetration 
of stellar radiation to photodissociate S0 2 , but not sufficient to ionize sulfur. Similarly, 
chemistry also affects the Fell and Sill lines, as there is substantial amounts of neutral 
iron and silicon in the denser regions (i.e. near the inner midplane) of the more massive 
disks. The SI line is optically thick for all the gas masses considered here and at higher gas 
masses, the line luminosity is approximately independent of mass. There is a slight increase 
in the line strength at 10~ 2 Mj as X-rays also contribute to the gas heating making the 
gas warmer, and a sharp decrease at lower gas masses as sulfur gets ionized due to lower 
column densities. The [SI]25.23/im line is stronger than the 56 /zm line by a factor of 3 — 10 
in all cases. The Sill and Fell line luminosities show a similar dependence with gas mass, 
becoming independent of mass the lines get optically thick for gas masses > 10~ 2 — 0.1 Mj. 
The emission from [Fel]24/zm, on the other hand, increases rapidly with gas mass (> 0.1 Mj) 
due to the increased abundance of warm dense Fel gas. The [Fel]24/zm line is the strongest 
line, with luminosities exceeding that of [SI]25.23;um for massive disks, making it an excellent 
tracer of gas in disks with masses > 0.1 Mj. However, its line strength depends on the gas 
phase abundance of Fe in these disk regions, which is uncertain. 

The [OI]63/im line is strong, and does not show much variation in strength with in- 
creasing gas mass, reflecting the fact that the [OI] line is nearly always optically thick. The 
emission is only weakly dependent on the dust mass, as it arises mostly from the surface 
regions of the disk where X-rays dominate the heating. The [OI]146/im line is weaker than 
the 63 jum line by approximately a factor of 10, but is strong enough to be detected in more 
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massive (and more dusty) disks. Emission from the [CII]158/xm line from r < 20 AU is 
too weak to be observed by Herschel or SOFIA, as the region of ionized carbon at the disk 
surface is too low in mass to be significant. CO emission increases with increasing gas mass. 
CO emission is optically thick for low J transitions and is dependent on the temperature and 
area of the emitting surface. The line strength increases for moderately high J transitions, 
whose wavelengths lie closer to the blackbody peak of the warm disk gas, and peaks at the 
optically thin 15 — >■ 14 transition at 187yum which arises mostly from within ~10 AU of the 
disk. 

Emission from various lines of H2O and OH is found to be strong for disks with gas 
masses > 0.1 Mj. In these disks, the densities are high and the abundance of these species is 
high making them strong coolants. OH is a strong coolant when the gas temperature is high 
(~few 100 K), and OH line emission mostly arises from the dense, hot inner disk where H 2 
photodissociates to OH more readily. Cooling by H 2 dominates in the dense equatorial 
regions of these massive disks at all radii, but the emission is spread in various lines in the 
mid-IR, only some of which lie in the the Spitzer Space Telescope band from 10—37 /im. The 
abundance of lines also implies that the cooling luminosity is distributed over many lines 
which reduces the strength of each individual line. Table 4 lists the strongest infrared line 
luminosities of the ortho and para-H 2 species and OH for a disk with a gas mass of 0.1 Mj 
and dust mass 10" 5 Mj, with all other parameters having their fiducial values. 



5.2. Variation of Spectra with Other Parameters. 

The position of the inner radial cutoff, is an important parameter which can alter the 
observed strength of emission lines. We have run cases with four different values of this inner 
radius equal to 0.3, 1.0, 3.0 and 10 AU from the central star, keeping all other parameters 
identical to our standard disk model. Disks with smaller inner holes have a small mass 
of hotter, denser gas which lies closer to the central star, and might be expected to have 
increased line luminosities as compared to disks with large inner holes. However, larger inner 
holes increase the gas mass heated and chemically altered by X-rays and UV, so the situation 
is in fact more complicated. The fluxes for all the lines (except Fel) are seen to drop when 
the inner radius is decreased from 1 AU to 0.3 AU. The higher disk densities in the inner 
disk cause sulfur to turn molecular, and reduce the ionization of iron and silicon. With lower 
abundances of S, Fe + and Si + and a smaller surface area of warm gas, their line luminosities 
correspondingly decrease. The line luminosity from the Fel line however increases due to 
a higher abundance of neutral iron. The strength of the H 2 lines is also seen to decrease 
slightly, and this is due to lower masses of > 100 K H 2 gas for the rj = 0.3 AU disk because 
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of effective self-shielding from UV and X-rays. When the disk inner radius is increased from 
3 AU to 10 AU, emission from all the lines except SI and H 2 S(0) increases again. A disk 
with such a large inner radius has a significantly low radial column density (~ 4 x 10 23 cm -2 
through the midplane) so that X-rays more effectively penetrate the disk, and now dominate 
gas-grain collisions as the main heating mechanism. The enhanced emission in most of the 
lines is due to a higher gas temperature due to X-ray heating through most of the disk. 
Emission from the [SI] lines is seen to decrease, however, and this is because of a higher 
ionization rate of sulfur which lowers the atomic sulfur abundance and lowers the effective 
area of warm [SI] gas relative to a disk with a smaller inner disk radius. 

Variations in the surface density power-law profile (Figure 13) have similar effects on 
the line intensities as changing ?v The amount of warm gas in the inner regions is sensitive 
to the power-law index a, and for a steeper power-law there is more mass in the inner 
region as compared to a shallower surface- density profile. We have considered three different 
power-law indices, a = 0.5, 1, and 1.5. For the largest value of a, there is a decrease in 
the strength of the [SiII]34.8/im and [FeII]26/xm lines, as the increased densities in the inner 
regions increase the shielding of these species from photoionization, and their abundances 
decrease. Correspondingly, the [Fel]24/xm line increases sharply for large a. 

The X-ray luminosity of the star was taken to be 10 -4 times the bolometric luminosity in 
the standard typical for a weak-lined T Tauri star or more evolved stars (Feigelson 

et al. 2003). However, this value may be higher or lower depending on the age of the star 
and other factors, and we have run cases where we assume an X-ray luminosity of 10~ 5 
to 10~ 3 times the bolometric luminosity. The X-ray spectrum is assumed to remain the 
same, peaking at around 2 keV and ranging from 0.5—10 keV. The X-ray luminosity has 
surprisingly little effect on the Spitzer Space Telescope lines as can be seen from Figure 14. 
The line luminosities seem largely insensitive to the amount of X-ray flux, and this can be 
understood by referring to the temperature structure (Figures 5 and 6) for the standard case. 
The X-rays only penetrate and heat a very small fraction of the mass of the disk, and in 
these regions, most of the heating is carried away as infrared emission from dust and by the 
[OI]63/im line. The strength of this line correlates with the X-ray flux as it comes mostly 
from the upper surface regions of the inner disk where the gas is warm and heated by X-rays 
(see Figure 6). 

The infrared line emission from gas is quite sensitive to the stellar flux, with stronger 
emission from gas disks around more luminous stars. We have run our disk models for two 
different stellar types, a G2 star and a K4 star. In both cases, all disk parameters were the 
same, and the X-ray luminosity in both cases was 10~ 4 the stellar luminosity. With lower 
X-rays and UV flux and with the dust grains around the less luminous K star being colder, 
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the gas temperature is much lower in the disk around the K4 star. Consequently, the line 
emission strengths are correspondingly lower (Figure 15). The [SI]25.23/im, [SiII]34.8/xm, 
[FeII]26/im and the [OI]63/im lines may still be strong enough for detection by the Spitzer 
Space Telescope and SOFIA for disks at d < 30 pc, especially if the X-ray luminosities are 
higher than average. 

Lowering the minimum grain size in the dust distribution increases line emission, as 
can be seen from Figure 16. In our standard case, the minimum grain size, a m j n was chosen 
to be 1 /xm, and we have run cases where the minimum grain size was 0.1 //m and 10/im. 
Smaller dust grains are less efficient at re-radiating stellar radiation, and get hotter as the 
grain size decreases. In addition, for a fixed dust mass and s = 3.5, the grain surface area is 
proportional to a~^ 2 , and therefore decreasing the surface area. Both of these 

effects raise the gas temperature, since the gas heating is largely by gas-grain collisions and 
this heating rate is proportional to (T^ust — T gas ) and the dust surface area. The resultant 
higher gas temperatures raise the total infrared line emission from the disk. 

In all our models, we did not allow the formation of H 2 on dust grains. It is not clear if 
hydrogen atoms stick on warm dust grains long enough to allow the formation of H 2 (Cazaux 
& Tielens 2002) and we have run the standard model where we allow the grain formation 
of H2. We find that there are no significant differences in the resultant spectrum of the disk 
in the two cases. Even when grain formation of H2 is set to zero, H2 is formed efficiently by 
reactions with H~ and by 3-body reactions. 



6. Discussion, Summary and Conclusions 

We have modeled the ~0.3 AU to 20 AU regions of "intermediate-aged" and "debris" 
disks with a range of dust masses 10 -7 — 10~ 4 Mj(or f D = L IR / L bo i = 10~ 4 — 10 _1 ) and 
gas masses 10~ 3 — 1 Mj. We define dust as grains with size a < 1 mm, and, assuming a 
size distribution n^ ust (a) oc a -3 ' 5 , we find that these disks are reasonably optically thin to 
starlight for M^ust < 10~ 4 Mj. Therefore, we have covered all of parameter space relevant to 
intermediate-aged and debris disks with sufficient gas and dust to be detected in IR emission 
lines in the next decade by observatories such as the Spitzer Space Telescope, SOFIA and 
Herschel. 

Although the dust is optically thin to stellar photons, the gas is optically thick at UV 
and X-ray wavelengths. Therefore, we find that the dominant gas heating mechanism for the 
bulk of the gas which is shielded from heating by UV or X-rays by gas in the inner regions 
or at the surface, is gas-dust collisions. In these regions the dust is heated by stellar photons 
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and colder gas is heated by collisions with warmer dust grains. However, at the surface of 
the disk and at the inner edge of the disk (assuming a hole or a gap between disk and star), 
the UV and X-rays heat the gas, and the gas temperature can exceed the temperature of 
even the smallest grains. We find that the gas temperature in the 0.3 AU to 20 AU regions 
of debris disks range from ~1000K to ~50 K. 

The chemistry in the innermost regions and at the surface of the disk is dominated by 
photoionization and photodissociation caused by UV and X-rays. However, in the cooler, 
shielded regions the chemistry is dominated, like interstellar cloud chemistry, by ion-neutral 
reactions with ionization by cosmic rays, UV and X-rays. We assume that at r < 20 AU, 
the grains are too warm {Td us t > 100 K) to support icy mantles or to allow surface catalysis 
of H 2 . However, H 2 is still efficiently formed by the reaction of H with H and by three- 
body reactions. In warmer regions {T gas > 300 K) neutral-neutral reactions with activation 
barriers become important, with the general result of enhanced production of OH and H 2 0. 
In addition, the reaction of O with H 2 to form OH can dominate the destruction of the 
H 2 molecule. Photodissociation of molecules by UV photons play a large role in molecular 
destruction and creation of ions such as Fe + and Si + . We find that sulfur is S + at the disk 
surface, S at intermediate gas opacities, and S0 2 in regions well shielded by overlying gas 
opacity. Because H 2 and CO so readily self-shield and are shielded by atomic C, we find that 
for cases with sufficient gas mass to be detected (M gas > 10~ 2 Mj) the bulk of hydrogen 
is H 2 and the bulk of carbon is CO in the disks. Unlike the case of interstellar clouds, the 
carbon becomes CO before the hydrogen becomes H 2 as one moves to greater depths from 
the disk surface, a result caused by the lack of production of H 2 on grain surfaces in the 
warm disk regions. 

Surprisingly, we find that the strongest mid-IR line is [SI]25.23//m. This results from 
the high value of the product of its abundance times its spontaneous decay rate (or A- value). 
We find that [SI] is a strong coolant in the inner disk for all our disk models with various 
input parameters. We have included detailed sulfur chemistry to determine the abundance of 
atomic sulfur in the disk, and emission is strong enough for detection even when substantial 
amounts of gas has turned molecular (mainly S0 2 ). The total mass of gas in disks is a little 
more difficult to obtain from measurements of the [SI]25.23yum line alone, as the amount of 
neutral atomic sulfur is affected by the chemistry of the disk. In addition, the [SI] line tends 
to be optically thick, and therefore not sensitive to the total gas mass. 

Other strong mid-IR lines include [SiII]34.8/xm, [FeII]26/im, [Fel]24/im and OH and 
H 2 rotational transitions. We assume depleted gas-phase silicon and iron abundances (see 
Table 5, Appendix A) similar to ice-free regions of the ISM. However, we have no iron chem- 
istry in our models other than ionization and recombination. We do include more chemistry 
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for silicon and consider some molecular species of Si, such as SiO and SiH 2 . However, the 
abundances of the silicon molecules are calculated to be very low, primarily because of pho- 
todissociation. In disks with low gas masses (M gas < 10~ 2 Mj), Si and Fe are ionized through 
most of the disk, and the line intensities we compute are probably reasonably accurate even 
though we do not treat the chemistry of these species in detail. In disks with higher gas 
masses, Si and Fe are shielded from photoionization and the [Fel]24/xm line becomes strong. 
The OH and H 2 lines become detectable at relatively high gas masses, M gas > 1CT 2 Mj. 

Pure rotational H 2 lines such as S(0) 28 /im and S(l) 17 /im are roughly an order of 
magnitude weaker than the [SI]25.23/xm line. The highest spectral resolving power of the 
Spitzer Space Telescope is R~600, and this modest resolution makes it difficult to detect weak 
mid-IR (10—37 /im) lines in the presence of strong(dust) mid-IR continuum. Detection of the 
relatively weak H 2 lines above the dust continuum requires a relatively large ratio (> 1000) 
of gas mass to dust mass. Such large ratios are possible at a transition epoch (~10 7 years) 
where much of the dust mass has coagulated into particles of sizes a>l mm, but the gas 
has not yet appreciably dispersed. 

Most of the Spitzer Space Telescope high spectral resolution emission arises from gas 
at r < 5 AU in our models, whereas longer wavelength (lower excitation) emission such as 
from low to mid-J CO lines and 01 63 /im has significant contribution from gas at r > 5 AU. 
The OI line arises from the entire disk and is a good probe of the conditions throughout the 
disk. While the mid-IR lines observable by the Spitzer Space Telescope are insensitive to the 
strength of the X-ray flux, the intensity of the [OI]63/im line does increase with the X-ray 
luminosity. SOFIA is capable of observing this transition, and may be able to detect smaller 
amounts of gas in disks around stars that are strong X-ray sources. Our models indicate 
that CO emission is strong enough to be detectable and it is a strong coolant in the outer 
cold regions of the disk. CO line luminosities are ~10~ 9 L in the 2 — >• 1 line, and ~10~ 7 
L in the higher 15 — ► 14 transition for our fiducial disk. Higher J transitions of CO in the 
far-infrared and the sub-millimeter wavelength regions can be easily detected by Herschel 
and can reveal the extent of the gas disk. 

Line ratios are also diagnostic of the gas characteristics. The [Fel]24/xm line is only 
strong when there is high shielding (or high columns and densities) in the inner disk regions. 
This requires either high gas masses, steep surface density distributions, small inner holes, 
or a combination of these. The Fel/SI ratio is especially suited to probing these parameters. 
Significant amounts of [OI]63/im and CO (J= 15 — >■ 14) 187 /im radiation emanate from 
regions r > 5 AU, in contrast to the SI or H 2 17 /im lines. Therefore, the ratios [OI]63/im/SI, 
[OI]63/im/H 2 17 /im, CO 187 /im/H 2 17 /im are sensitive to the relative masses of gas at 
r < 5 AU compared to that at r > 5 AU. In other words, these ratios are sensitive to the 
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steepness of the surface density distribution. 

The overall strength of the mid-IR lines primarily depends on the mass of warm (T > 100 
K) gas, which typically lies within 5 AU of the central star. Disks with more total gas mass, 
with smaller inner holes, or with steeper surface density distributions, produce stronger mid- 
IR emission. Because gas-dust collisions mainly heat the bulk of the gas, the line emission also 
increases with either greater dust mass, or, holding dust mass fixed, with smaller minimum 
grain size a m j n . In the latter case, the smaller size grains are warmer and have greater surface 
area than their larger counterparts, and therefore heat the gas more vigorously. 

As little as 10~ 3 Mj of warm gas at r < 5 AU can be detected by the Spitzer Space 
Telescope in the [SI]25/xm line for sources at ~30 pc from the Sun. However, to better 
diagnose the physical characteristics of the gas requires the detection of several lines, and 
this requires roughly an order of magnitude greater gas mass (~10~ 2 Mj) at this distance. 
We emphasize that the line luminosities in the strongest lines ([SI], [Sill], [Fell] and [01]) 
increase very slowly or even decrease with increasing gas mass (see Fig. 10). This is because 
the gas temperature in these disks is largely controlled by the dust temperature, which is 
independent of gas mass, and the gas emission lines at high gas mass become optically thick 
so that their luminosities saturate. Our results show that at the distance of Taurus (140 pc), 
even 1 Mj of gas in intermediate- aged optically thin dust disks could be difficult to detect. 
However, typical disks in Taurus are more likely to be younger and therefore abundant in 
small grains (a min ~ 50A). A younger central star is also more luminous in FUV. This will 
lead to an increase in both collisional gas-grain heating (due to the larger surface area of 
small grains) and grain photoelectric heating which will make the gas warmer and increase 
gas line emission significantly. 

Assuming that we know the distance to the source, the X-ray and FUV luminosities of 
the central star, and the dust mass and dust distribution from the IR and sub-mm/millimeter 
continuum observations, we can readily constrain the gas properties by comparing the ob- 
served fluxes in various emission lines with our models. Ultimately, if the lines are spectrally 
resolved by high resolution instruments on SOFIA or Herschel, the line profiles will also help 
constrain the gas dynamics and spatial distribution. 

To conclude, we find exciting prospects for the detection by the Spitzer Space Telescope, 
SOFIA or Herschel of infrared line emission from the gas in 0.3 AU to 20 AU regions of 
intermediate-aged disks at distances < 150 pc. The core accretion model of giant planet 
formation strongly suggests an epoch (t < 10 6 5 — 10 7 years) when the dust in the 1 AU— 
20AU region is relatively optically thin, due to the growth of several earth-mass cores, but 
where > 1 Mj of gas persists for an extended period. Therefore, the observations of line 
emission coupled with the models presented here will provide a measure of the evolution of 
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gas in disks as they transition from planet-forming disks to the older, gas-free, debris disks 
and a test of the core-accretion scenario for gas giant planet formation. The presence of even 
small amounts of gas at times t ~ 10 7 to 10 8 years has important effects on planet formation 
in the 0.3-20 AU region. 
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Table 1: Input Parameters for Standard Disk Model 



Gas mass 


KT 2 Mj 


Dust mass 


10^ 5 Mj 


Inner radius r« 


1 AU 


Minimum grain size a m i n 


1.0 /im 


Outer radius r Q 


20 AU 


Maximum grain size a max 


1.0 mm 


Surface density power law index a 


1.0 


Grain size dist. index s 


3.5 


Stellar age (for FUV flux) 


10 7 years 


Grain composition 


Am. Olivine 


Stellar X ray Luminosity 


4.0 10 29 erg s" 1 


Stellar bol. luminosity 


4.0 10 33 erg s" 1 


Stellar temperature T* 


5780 K 







Table 2: Infrared Emission Lines, their Wavelengths and Observing Instruments. 



Species A(/im) Transition Telescope instruments 



H 2 S(2) 


12.3 






Spitzer 


IRS, SOFIA EXES 


H 2 S(l) 


17.0 






Spitzer 


IRS, SOFIA EXES 


H 2 S(0) 


28.2 






Spitzer 


IRS, SOFIA EXES 


OH 


15-200 






Spitzer 


IRS,SOFIA EXES, GREAT, SAFIRE, Herschel SPIRE 


H 2 


15-200 






Spitzer 


IRS,SOFIA EXES, GREAT, SAFIRE, Herschel SPIRE 


Fel 


24.0 


5 D 3 


- 5 D 4 


Spitzer 


IRS, SOFIA EXES 


SI 


25.2 


3 Pi 


- 3 P 2 


Spitzer 


IRS, SOFIA EXES 


Fell 


26.0 


6 D 7/2 


- 6 D 9/2 


Spitzer 


IRS, SOFIA EXES 


Sill 


34.8 


2 P 3 /2 


- 2 Pi/ 2 


Spitzer 


IRS 


Fel 


34.2 


5 D 2 


- 5 D 3 


Spitzer 


IRS 


Fell 


35.4 


6 D 5/2 


- 6 D 7/2 


Spitzer 


IRS 


SI 


56.6 


3 Po 


- 3 Pi 


SOFIA 


FIFI-LS 


OI 


63.1 


3 Pi 


- 3 P 2 


SOFIA 


FIFI-LS, GREAT, Herschel PACS 


OI 


145.6 


3 Po 


- 3 Pi 


SOFIA 


SAFIRE, GREAT, Herschel SPIRE 


CII 


157.7 


2 P 3/2 


- 2 Pi/ 2 


SOFIA 


SAFIRE, GREAT, Herschel SPIRE 


CI 


369.0 


3 P 2 


- 3 Pi 


SOFIA 


CASIMIR, SAFIRE, Herschel HIFI 


CI 


609.2 


3 Pi 


- 3 Po 


SOFIA 


SAFIRE, Herschel HIFI 



Table 3: Parameter Range for Disk Model 



Gas mass 


10" 3 - 1 Mj 


Dust mass 


10" 7 - 10" 4 Mj 


Inner radius 


0.3-10 AU 


Minimum grain size a m i n 


0.1-10.0 /im 


Outer radius r Q 


20 AU 


Maximum grain size a max 


1.0 mm 


Surface density power law index a 


0.5-2.0 


Grain size dist. index s 


3.5 


Stellar X ray Luminosity 


10-2 _ iQ-5 Lbol 


Stellar temperature T* 


4500K and 5780 K 
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Table 4: Line luminosities of OH and H2O lines for a disk with a gas mass of 0.1 Mj and 
dust mass 10~ 3 Mj, with all other parameters at their fiducial values. 



Species 


A(/im) 


Line lum. (in L Q ) 


Species 


A(/im) 


Line lum. (in L Q ) 


OH 


30.54 


9.8xl0~ 8 


OH 


30.63 


l.OxlO- 7 


OH 


33.80 


1.4xl0~ 7 


OH 


33.93 


1.4xl0- 7 


OH 


37.01 


1.8xl0~ 7 


OH 


37.08 


1.8xl0~ 7 


OH 


38.84 


7.7xl0~ 7 


OH 


38.88 


6.5xl0^ 7 


o-H 2 


21.15 


2.6xl0~ 7 


P-H 2 


21.17 


2.6xl0~ 7 


P-H 2 


21.18 


3.0xl0~ 7 


P-H 2 


21.83 


2.8xl0~ 7 


o-H 2 


21.85 


3.1xl0~ 7 


P-H 2 


21.89 


2.7xl0~ 7 


P-H 2 


23.19 


3.0xl0- 7 


o-H 2 


23.90 


3.2xl0~ 7 


o-H 2 


25.17 


2.9xl0~ 7 


P-H 2 


25.22 


3.1xl0- 7 


o-H 2 


25.94 


4.0xl0- 7 


P-H 2 


25.98 


3.1xl0- 7 


P-H 2 


27.03 


3.1xl0- 7 


o-H 2 


27.99 


3.1xl0- 7 


o-H 2 


28.59 


2.7xl0~ 7 


P-H 2 


28.59 


2.7xl0- 7 


P-H 2 


28.91 


2.3xl0~ 7 


P-H 2 


29.36 


2.7xl0~ 7 


o-H 2 


29.84 


3.0xl0- 7 


P-H 2 


29.88 


2.9xl0~ 7 


P-H 2 


30.47 


2.6xl0~ 7 


P-H 2 


30.53 


2.7xl0~ 7 


o-H 2 


30.53 


2.8xl0~ 7 


o-H 2 


30.87 


2.5X10" 7 


o-H 2 


30.90 


3.4xl0~ 7 


P-H 2 


31.74 


2.5xl0~ 7 


P-H 2 


33.00 


2.9xl0~ 7 


o-H 2 


33.01 


3.0xl0~ 7 


P-H 2 


34.40 


2.0xl0- 7 


o-H 2 


34.55 


2.7xl0~ 7 


o-H 2 


35.43 


2.7xl0~ 7 


P-H 2 


35.47 


3.2xl0~ 7 


o-H 2 


35.67 


2.2xl0~ 7 


P-H 2 


35.90 


2.9xl0" 7 


P-H 2 


36.21 


2.7xl0~ 7 


o-H 2 


66.44 


2.6xl0~ 7 


o-H 2 


75.38 


3.6xl0" 7 


P-H 2 


100.98 


8.1xl0~ 7 


o-H 2 


108.07 


3.3xl0~ 7 


P-H 2 


138.53 


5.9xl0~ 7 


o-H 2 


179.53 


2.2xl0~ 7 
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A. Chemistry 

We adopt gas phase abundances from Savage & Sembach (1996) towards the line of sight 
through a diffuse cloud to the star ( Oph. The abundances are listed in Table 5. Table 6 lists 
the atomic, ionic and molecular species included in our chemical code. Reaction rates are 
mostly from the UMIST'99 database, except for some reactions which are listed in Table 7. 
Apart from those listed here, X-ray photoionization rates are calculated as described in 
§D.3, and reactions involving PAHs and vibrationally-excited H 2 are from TH85. The UV 
photoionization rates of atoms and molecules are calculated using their photoabsorption 
cross sections multiplied by the stellar flux and integrated over frequency. 



B. X-ray Heating and Chemistry 

B.l. X-ray Spectrum 

The stellar X-ray spectrum is assumed to be similar in form to the wTTS in Figure 2 
of Feigelson & Montmerle (1999). A broken power law is fit to the data for the counts G(in 
counts s _1 keV -1 ) within the energy range 0.5-10 keV, which then gives the luminosity in 
erg s" 1 keV" 1 , L(E) = G(E)E. 

L(E) = 1.2L X E- 175 (E > 2keV) (Bl) 
= 0.18L X E(E < 2keV) (B2) 

(B3) 

where Lx is the total X-ray luminosity of the star in erg s _1 , and the energy E is expressed 
in keV. The X-ray energy flux is given by F(E) = L(E)/4:iiD 2 , where D = (r 2 + z 2 ) 1 ^ 2 is 
the distance from a point in the disk to the star. The heating due to X-rays is 

/•10 

H x = F(E)e- Tx a(E)n H f h dE erg cm" 3 s" 1 (B4) 

io.5 

where cr(E) is the total gas X-ray photoabsorption cross section per H nucleus, Tx is the 
X-ray optical depth given by a(E)N H , N H is the total column density of the gas towards 
the source, is the number density and fh is the fraction of absorbed energy that goes into 
heating the gas. We also include secondary heating of X-rays due to photoionization from 
the n = 2 state of atomic hydrogen and follow exactly the treatment of Shang et al.(2002). 
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B.2. X-ray Photoabsorption Cross Section a(E) 

We use the model of Wilms et al. (2000) for the X-ray photoabsorption cross section per 
H nucleus. They obtain the total photoionization cross section of the ISM by summing the 
cross sections of the gas, molecules and grains. They use cosmic abundances for the various 
elements, and add the contribution from grains. Our grains are large and typical X-ray 
penetration depths are ~0.1//m, and therefore X-rays do not penetrate even our smallest 
dust grains. We use only the gas-phase abundances in our calculation of the cross-section, 
and we adopt abundances from Savage & Sembach (1996) through a diffuse cloud towards 
the star ( Oph. We correct Figure 1 of Wilms et al.(2000) to correspond to our depleted 
abundances and then fit a power law to find the absorption cross section. We obtain 

/ F \ -2-594 

a(E) = ,. 2 x 10 -» (_) (B5) 
Compare to the fit to the Wilms et al.(2000) cross-section 

/ rp \ -2.485 

a(E) = 2.27 x 10~ 22 (B6) 



B.3. X-ray Photoionization Rates of Atoms and Molecules 

The absorption cross section of a species i for ionization is approximated to be of the 
functional form 

(E \ a 
IkevJ 

We use the partial photoionization cross sections of Verner & Yakovlev (1995) and the fitting 
formulae they provide (Table 1 of their paper). These formulae for the cross sections in the 
range 0.5-10 keV are then fit again with a power law of the above form for each neutral 
atomic species. In some cases, the cross section is fit with a broken power law within this 
range. Table 8 lists the values of a 0i in barns, a and the energy ranges in keV where the fits 
are valid. The primary rate of photoionization is given by the integral of the photon number 
flux and the cross section 



d = 6.25 x 10 8 / ai (E)l-^e~ Tx( ^dE 

Jo.5 E 



where F(E) is in units of erg s" 1 keV -1 , and the energy E is in keV. The attenuation factor 
Tx has the same definition as before and is equal to the product of the total photoabsorption 
cross section and the column density towards the source. 
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Secondary ionizations are important for hydrogen and helium. The fraction of primary 
electron energy going into secondary ionization is (Shull & van Steenberg 1985), 

13.6 - 41 ^ L76 



$ H = 0.391 ( 1 - — ) « 0.35 



24 g0.46\ L66 



<f> He = 0.0554 ^1 - — j « 0.05 

where E is now in eV. The number of secondary photoionizations produced by a primary 
electron of energy E is equal to -E$#/13.6 and -E , $# e /24.6 for hydrogen and helium. The 
number of secondary ionizations by a photon of energy E is therefore 25.7(£'/lkeV) for 
hydrogen and 2.0(£'/lkeV) for helium. The total photoionization rate for hydrogen and 
helium are therefore, 

( TH = 6.25 x 10 8 / a i {E) 1 -^-e- Tx{E \l + 25.7E)dE s" 1 
io.5 E 

( THe = 6.25 x 10 8 / o i {E) 1 -^-e- Tx{ - E \l + 2.0E)dE s" 1 
Jo.5 E 

where E is in units of keV. We ignore secondary ionizations of all other species. Photoion- 
ization rates of molecules are assumed to be the sum of the rates of their constituent atoms, 
e.g., (oh — Co + (h e tc. These integrals are numerically evaluated for each species with the 
expressions for the cross sections and the X-ray flux, as a function of column Nh towards 
the star. 



C. Gas Heating Mechanisms 

Gas-grain collisions Dust grains in the optically thin disk are heated by the stellar radi- 
ation field and collisional gas-grain energy transfer is often the dominant heating mechanism 
for the gas. The dust grain temperature is size-dependent, smaller grains are less efficient at 
re-radiating stellar light and are thus warmer than larger grains. Collisions of gas molecules 
with dust grains can either heat or cool the gas, depending on the gas and dust grain tem- 
peratures. Our treatment is similar to that of Hollenbach & McKee (1979, hereafter HM79) 
but adapted for an arbitrary grain size distribution. The net heating (or cooling) due to 
collisions is given by 

Q>max 

r\ = a T n H Vgas ^2 na2n dust(a)C2kBT dust (a) - 2k B T gas ) ergs cm" 3 s" 1 (CI) 
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where cut ~ 0.3 is the thermal accommodation coefficient (Burke & Hollenbach 1983), uh 
the gas hydrogen nucleus number density, v gas the thermal velocity of the gas, a the dust 
grain radius, that ranges from a min to a max , nd ust (a) the number density of particles of size a, 
and T dust (a) and T gas are the temperatures of a dust grain of size a and the gas respectively. 

FUV pumping of H 2 and collisional de-excitation H 2 molecules absorb FUV photons 
from the star and from the ISRF and are excited to a higher vibrational state. Subsequent 
de-excitation by collisions with H atoms and ground state H2 molecules deposits the excess 
energy as heat into the gas. H 2 molecules can also be directly photodissociated from the vi- 
brationally excited state, which again leads to heating. (See TH85 for a detailed discussion). 
The heating due to FUV pumping and collisional de-excitation is given by 



where n(H 2 ) and ^(rD;) are the number densities of H 2 and vibrationally excited H 2 respec- 
tively, Pi(t) and /3 2 (r) the corresponding self-shielding factors and Gn.26 is the local photon 
energy flux in the range 11.26 — 13.6eV (from the star and the ISRF, attenuated by dust 
and other gas species), normalized to the Habing flux in the same band (2.26 x 10~ 4 erg 
cm -2 s _1 ). We use self-shielding factors for CO from van Dishoeck & Black (1988) and for 
H 2 from Draine k Bertoldi (1996). 

X-ray heating X-rays incident on gas photoionize atoms and molecules, and the primary 
electrons thereby generated are energetic enough to cause multiple secondary ionizations. 
The heating term arises from the thermalization of some of the energy of these secondary 
electrons, and is proportional to the ionization rate due to X-rays. The heating due to X-rays 
is given by 



where F(E) is the X-ray photon flux, Nh is the hydrogen nucleus column density of gas 
towards the source, cr(E) is the X-ray absorption cross section of the gas per hydrogen 
nucleus (see Appendix B), and fh is the fraction of energy that is converted to heat, equal 
to 0.1 for atomic gas and 0.4 for molecular gas. (See Maloney et al. 1996 for details.) We 
also include indirect heating by X-rays due to atomic hydrogen ionization, the recombining 
H + resulting in a substantial n = 2 population of H because of the trapping of Lyman a 



r 2 = 



3.18 x 10- 12 n(H;)(10- 12 T a5 exp(-1000/T)n(/J) + 
1.4 x io~ 12 T a5 exp(-18100/(T + 1200))ra(H 2 )) + 
G n . 26 (1.36 x 10- 23 (n(H 2 )/3 1 (r)n(H;))/3 2 (r)) + 
5.8 x 10" 24 n(H;)) ergs cm" 3 s" 1 



(C2) 




(C3) 
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photons, followed by photoionization or collisional de-excitation of n = 2 hydrogen atoms. 
(See Shang et al. 2002.) We ignore X-ray scattering from the surface layer into the disk, as 
we estimate the contribution from this process to be negligible. 

Grain photoelectric heating Energetic electrons emitted from dust grains due to the 
absorption of stellar FUV photons can collide with and heat the gas (Bakes & Tielens 1994). 
The efficiency of this process is dependent on the size of the dust grain and decreases with 
increasing grain size. As we consider large dust grains in our disk models, this process 
is not as efficient as in the ISM where small submicron sized dust grains are abundant. 
The grain photoelectric heating term has been revised to include the effects of a cooler 
radiation field (stars of lower spectral type) and increased grain sizes. We use corrected 
grain photoelectric efficiencies from Weingartner & Draine (2001) and Weingartner (2002, 
private communication) for our calculations of grain photoelectric heating in disks. 



where e is the photoelectric heating efficiency and is a function of the grain size a, the equiv- 
alent blackbody temperature of the stellar FUV flux T FUV) and G and n e (which determine 
the grain charge). The FUV flux spectrum from the star (Figure 1) has a characteristic 
temperature Tpjjy ~ 10, 000 K. Go is the local stellar flux between 6.0—13.6 eV, normalized 
to the Habing flux. A similar term for the ISRF is also included, where Tpuv is replaced by 
the temperature of the ISRF, 30000 K, although this contribution is generally negligible. 

Drift heating Dust grains and gas molecules are subject to different forces in the disks, 
with radiation pressure acting only on dust and gas being constrained by thermal pressure 
gradients. They could have different disk velocities and orbits, and this causes a viscous 
heating term which could be important in dense disks (Kamp & van Zadelhoff 2001, Takeuchi 
and Artymowicz 2002). We add drift heating to our disk models, though we do not solve for 
the disk orbital dynamics of the dust or gas. The heating due to drift is given by 



where Vdrift(a) is the drift velocity between gas and a dust particle of size a. We use the 
equations of Takeuchi and Artymowicz (2002) for the calculation of the drift velocity. 




(C4) 




(C5) 
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Heating due to formation of H 2 We include heating due to the formation of H 2 in 
our model. The formation of H 2 either on grain surfaces or by reactions of H atoms with 
H~ leads to heating of the gas. The grain formation mechanism provides ~0.2 eV of kinetic 
energy going into the newly formed H 2 molecule, and ~4.2 eV of energy initially deposited 
as rovibrational internal energy of the molecule (Hollenbach, Werner & Salpeter 1971). The 
reaction with H~ is an exothermic reaction which liberates ~3.53 eV of energy and we 
assume that most of this energy is initially in rovibrational internal energy of the molecule. 
Generally in disks, the gas densities are much greater than the critical density n cr to de- 
excite the vibrationally excited H 2 molecule. Therefore, the rovibrational energy is quickly 
converted into heat. Grain formation of H 2 is not included in our standard runs, and this 
heating mechanism is only included for the case where we allow formation of H 2 on grains. 
Heating due to formation of H 2 by H is given by 

T 6 = 7.34 x 10- 21 n H 2 <H-)x{H) erg§ cm _ 3 s _, 
1 + n cr /n H 

where 

iq6 T -1/2 

n _ 1U 1 cm -3 cq 7 \ 

cr 1.6x(#)exp[-(400/T) 2 ] + 1.4x(# 2 ) exp[-12000/(T + 1200)] V ; 

and x(H),x(H~), and x(H 2 ) are the abundances of H, H~ and H 2 respectively (HM79). 
Heating due to formation of H 2 on dust grains is given by 

T 7 = 4.8 x 10- 29 n H 2 f dust x(H)(0.2 + - —. ) ergs cm" 3 s" 1 (C8) 

1 + n cr /n H 

and fdust is a factor that accounts for the change in dust surface area per H nucleus caused 
by changes in the size distribution and the gas-to-dust mass ratio as compared to interstellar 
values. 



Cosmic ray heating and other heating mechanisms Apart from the processes dis- 
cussed above, we consider cosmic ray heating (r 8 ), photoionization of neutral carbon (r 9 , see 
TH85), and heating by the excitation of [OI]63;um upper level by the IR continuum radiation 
field followed by the collisional de-excitation of this level (ri , see Appendix D). 

T 8 = 1.5 x lO- n (n H (x(H) + 2x(H 2 )) ergs cm" 3 s" 1 (C9) 

where ( is the primary cosmic ray ionization rate. 

T 9 = 4.16 x 10- 29 n H x(C)G n .26 exp(-3A v - Tn. 26 )/3 C ergs cm" 3 s" 1 (CIO) 



-42- 



where x(C) is the abundance of neutral carbon, rn. 26 is the optical depth due to gas absorp- 
tion towards the star in the FUV band that photoionizes carbon, and 0c is a self-shielding 
factor (van Dishoeck & Black 1988). To this term we also add similar terms for the contri- 
bution from the ISRF, and use optical depths and self-shielding factors due to gas columns 
vertically up and down from the spatial position in the disk. 

D. Gas Cooling Processes 

Gas cooling occurs by gas-grain collisions when the gas is warmer than the dust (Eq. CI) 
and by collisional excitation followed by line radiation of gaseous atoms, ions and molecules. 
We essentially follow the treatment of TH85, in which the cooling due to a transition from 
level i to level j is given by 

A i:j = riiAijhvijPescfaj) (1 - P(vij)IS{yij)) , (Dl) 

where rii is the population density of the species in level i, Aij is the spontaneous transition 
probability, hu^ is the energy of the emitted photon, is the optical depth in the line and 
Pesd T ij) is the associated escape probability. P(vij) is the background radiation term which 
consists of the sum of the cosmic microwave background radiation and the infrared emission 
from the dust grains of different temperatures at the wavelength of the transition. The 
source function S{fij) is determined by solving for the level populations assuming statistical 
equilibrium. We generally assume that radiation can escape from any point in the disk either 
vertically up or down or radially out the inner radius. We calculate the column densities 
in each level of the cooling species in all three directions. From the column densities we 
determine the optical depth in the line and the escape probability of a photon, which then 
gives the cooling rate in a given line. If the gas temperature is lower than that of dust, 
and the background radiation stronger than the source function, the transition can heat 
the gas as the term 1 — P(uij)/S(vij) becomes negative. In this case, a background photon 
is absorbed by a particular species and subsequent collisional de-excitation transfers the 
excitation energy to gas heating. 

We treat CO and H 2 cooling in a detailed multi- level calculation but treat OH and H 2 
somewhat differently. Analytical approximations to multi-level line cooling calculations of 
OH and H 2 are derived for use in the gas disk code. The full level population calculation 
was prohibitively expensive to include in code and instead, we derived analytic fits for the 
total cooling similar to HM79. For H 2 0, we used collisional rate data for approximately 170 
levels for both the ortho and para species, and new parameters (refer HM79) that best fit 
the data, o = 2.0 x 10~ 15 cm- 2 , E /K=23.0 K and A = 6.0 x 10~ 3 s _1 . The collisional 
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rate data for OH for the first 45 levels are used to obtain the values, a = 8.0 x 10 16 cm 2 , 
E /K=5A K and A = 7.6 x 10" 4 s~\ 
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Table 5: Elemental gas phase abundances (from Savage & Sembach 1996) 



Element 


Abundance 


He 


0.1 


C 


1.4 xl(T 4 





3.2 xl(T 4 


Si 


1.7 xl(T 6 


Mg 


1.1 xl(T 6 


Fe 


1.7 xl(T 7 


S 


2.8 xl(T 5 



Table 6: Chemical species included in the models 



H 


He 


c 





H 2 


o 2 


OH 


CO 


H 2 


H+ 


He+ 


c+ 


0+ 


OH+ 


CO+ 


H 2 0+ 


HCO+ 


H 3 0+ 


H 2 + 


H 3 + 


CH+ 


CH 2 + 


CH 3 + 


CH 


CH 2 


Mg 


Mg+ 


Si 


Si+ 


SiH 2 + 


SiH 


SiO 


Fe 


Fe+ 


S 


S+ 


SiO+ 


HOSi 


H 




CH 4 


CH 3 CH 4 + 


HS+ 


H 2 S + 


H 3 S+ 


HS 


H 2 S 


cs 


CS+ 


HCS+ HCS 


OCS 


OCS+ 


HOCS+ 




SO 


SO+ 


HSO+ 


so 2 


S0 2 + HS0 2 + 


H 2 CS 


H 2 CS+ 


H 3 CS+ 




s 2 + 


S 2 H+ 

















Table 7: Reaction rates where they differ from the UMIST'99 database. 

H 3 0+ + e" TTo ' il 1.0 x 10 (i V • 2 

H 3 + + e" -> H 2 + H + O 1.73 x lO^T' 1 / 2 

H 3 + + e" -> H 2 + OH 3.3 x m^T' 1 / 2 

H 3 0+ + e" -> OH + H + H 9.5 x lO^T" 1 / 2 

H 3 + + e" -> H 2 + H 2.9 x W^T' 1 / 2 

H 3 + + e" -> H + H + H 9.0 x lO^T' 1 / 2 

H + H + H -> H 2 + H 8.8 x 10" 33 

H 2 + H + H -> H 2 + H 2 2.8 x iQ- 31 T- a6 
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Table 8: Fits to photoionization cross sections where the cross section is given by <Ji(E) = 
a 0i {E/lkeV) a and E x < E < E 2 . 



Atom 


<7oi barns 


a 


E x keV 


E 2 keV 


H 


11.55 


-3.4 


0.5 


10 


He 


43.12 


-3.3 


0.5 


10 


C 


4.56 xlO 4 


-3.01 


0.5 


10 





7.34 xlO 3 


-3.125 


0.5 


0.537 




1.28 xlO 5 


-2.9 


0.537 


10 


Mg 


5.93 xlO 4 


-2.38 


0.5 


1.312 




5.5 xlO 5 


-2.76 


1.312 


10 


Si 


1.00 xlO 5 


-2.44 
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Fig. 1. — The adopted stellar spectrum for a 10 Myr old G2 star, consisting of an age- 
dependent UV flux added to a Kurucz model for the same spectral-type. The dashed line is 
a spectrum for a K star that includes chromospheric emission. 
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Fig. 2. — Disk surface for the standard case at a density equal to 10 the density in the 
midplane, and contours that indicate surfaces where the vertical gas column densities to the 
surface are 10 21 and 10 22 cm -2 . 
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Fig. 3. — Heating terms as a function of radius in the disk midplane for the standard 
disk model with M gas = Mj and M dust = 10~ 5 Mj. F 1 is the heating due to gas-grain 
collisions, T 3 is X-ray heating, T 4 grain photoelectric heating, T 6 heating due to H 2 formation 
and T 8 cosmic ray heating. See Appendix C for details. 
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Fig. 4. — Cooling due to different chemical species as a function of radius in the disk midplane 
for the standard disk model. 
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Fig. 5. — Gas and dust temperatures in the midplane of the disk are shown as (labelled) 
solid lines, whereas the dashed line shows the gas number density in the midplane for the 
standard disk model. 
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Fig. 6. — Gas temperature (solid line) and density (dashed line) as a function of z at a radial 
position of 2 AU and the temperature at 10 AU for the standard disk model. The dominant 
heating agents and the coolants in each region are marked on the temperature curve. 
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Fig. 7. — Abundances of important species as a function of z at a radius of 2 AU for the 
standard disk model. All abundances, except for OH and H 2 0, have been normalized to 
their maximum possible values. 



-53- 



1.2 



1 - 



2 0.8 

O 



0.6 



0.4 



i 1 1 r 



SI 



25 



Fell 



1 

0.5 


-0.5 



~i i r 



i i r 




J I I I I L_ 



30 



40 



-0.3 
-0.4 



LH 2 's(2) 


iii- 


: , 1 , 


1 l~~-4_ 



12 



12.5 




H 2° Sill 



H 2 S(0) 



30 



H 2 H 2 



35 



17.5- 




40 



Fig. 8. — Mid-infrared spectrum in the 24—40 /xm wavelength region showing the dust 
continuum and dominant gas emission lines for the standard disk model, assuming a distance 
to the disk of 30 pc and a spectral resolving power (A/AA)=600, typical of high resolution 
Spitzer observations. The dashed line indicates a 5% line-to-continuum ratio above which 
lines may possibly be detected. The insets show the complete spectrum from 10 — 40/im, 
and the S(l) and S(2) lines of H 2 . 
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Fig. 9. — Mid-infrared spectrum (as in Figure 8), for a disk with an inner radius of 10 AU, 
and with all other parameters the same as the standard model. 
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Fig. 10. — Mid-IR line luminosities as a function of disk gas mass, for different dust masses 
(10~ 5 Mj upper line, 10 -6 Mj middle line, and 10~ 7 Mj lower line) in the disk. The dashed 
lines show the luminosities of the [SI]56/xm line and the [Fel]24/zm line. 
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Fig. 11. — Detectability plots for three representative gas emission lines H 2 S(l), [SI]25.23/xm 
and [Sill] 34.8/xm, by the Spitzer Space Telescope (IRS instrument with R=600) for a stan- 
dard model disk at 30 pc. The lines show the minimum detectable gas mass for different 
dust masses in the disk for each line. For M dust < 1CT 5 Mj, detection is flux limited (3 x 1CT 8 
L for A ~ 25/im). For M dust > 10~ 5 Mj, the line-to-continuum ratio is the limiting factor. 
We assume that a line-to-continuum ratio of 5% is sufficient for detection. 
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Fig. 12. — Gas line luminosities for disk models with differing inner radii, Tj, and other 
parameters as in the standard case. 
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Fig. 13. — Changes in the line luminosities due to variations in the slope of the power law 
index a of the surface density distribution. 
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Fig. 14. — Gas line luminosities for disks around a G star with different X-ray luminosity. 
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Fig. 15. — Gas line luminosities for disks with parameters as in the standard disk model but 
around a G star(T* = 6000) and a K star(T* = 4500). 



- 61 - 





-o 




i i i i 
SI 25.2/u.m 

1 1 35.4/xm 


i i i i i i i 


i I i 


i 

— 




-6 




Fell 26/xm 


______ 












H 2 S(l) 






- 


o 

J 


-7 




H 2 S(2) 

H 2 S(0) 








o 


-8 




i i i i 


i i i i i i i 


ili 


- - 

iil~ 




-1.5 


-1 


-0.5 

L °g a mm (^m) 


0.5 


1 




-5 
















i i i i 


I I I I I I I 


i I i 


i 




-6 




01 63 /_m 

SI 56 fim 








o 


-7 
-8 




01 146 /xm 








Log L 




CO 187 fim 

i i i lii 


l l l l l l l 


ili 


ill - 




-1.5 


-1 


-0.5 


0.5 


1 



Log a min (jj.m) 

Fig. 16. — Gas line luminosities for disks with dust distributions where the minimum grain 
size a m j„ is varied. 



